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(54) Wheel information estimating apparatus 

(57) A wheel information estimating apparatus, 
which includes: a detecting device (10, 12, 18, 56, 140) 
for detecting a wheel motion state relating to a movement 
of a tired wheel of a motor vehicle; an estimating device 
(52, 53, 82, 84, 86, 94. 95. 160. 1 62) for estimating wheel 
information relating to the tired wheel, on the basis of the 
wheel motion state detected by the detecting device, and 
a modifying device (64, 80, 92. 60, 164) for modifying an 
element relating to estimation of the wheel information, 
so as to improve an accuracy of estimating the wheel 
information. The element to be modified is selected from 
a movement of a center of rotation of the tired wheel as 
viewed in a plane in which the wheel is rotated, a signal 
applied from the detecting device to the estimating 
means, and an internal parameter used by the estimating 
means for estimating the wheel information. 
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Description 

B C KG ROUND OF THE INVENTION 

5 F d of the Invention 

The present invention relates to an apparatus adapted to estimate a certain kind of information relating to a tired 
wheel (hereinafter referred to as "vehicle wheel" or "wheel") of a motor vehicle, on the basis of a parameter indicative 
of the motion of the wheel, and particularly to techniques for improving the accuracy of estimating the wheel information. 

10 

Discussion of Related Art 

For controlling a motor vehicle, it is useful and significant to known the motions and changing characteristics of the 
wheels of the vehicle, because the wheels are the only components at which the vehicle contacts a road surface. How- 

rs ever, it is not easy to directly detect the motions and the changing characteristics of the vehicle wheels during running 
of the vehicle, since the wheels are attached to the body of the vehicle through a suspension system such that the 
wheels are displaceable relative to the vehicle body, and since the wheels are rotated. 

For a motor vehicle equipped with wheels having tires that are inflated by compressed air so as to maintain a suitable 
degree of elasticity, for example, the air pressure in the tires has an important effect or influence on the operating char- 
ge acteristics of the wheels. Accordingly, it is desirable to detect the air pressure levels of the wheel tires, but the detection 
of the tire pressure is also difficult. 

The air pressure of a wheel tire may be detected by a pneumatic pressure sensor attached thereto. However, it is 
not easy to transmit the output of the sensor on the rotating tire to a desired position on the vehicle body while the vehicle 
is running. It is also possible to determine the tire pressure on the basis of a change in the distance between the road 

25 surface and the vehicle body, which change occurs due to deformation of the wheel tire as the tire pressure is lowered. 
Further, it is possible to estimate a decrease of the tire pressure of one or more wheels as compared with the tire pressure 
of the other wheels, depending upon the detected rotating speeds of the wheels, based on a fact that the tire pressure 
decreases with an increase in the rotating speed. However, these methods do not assure accurate detection of the tire 
pressure. In particular, the latter method does not permit the detection of a decrease of the tire pressure of any wheel 

30 if the tire pressure levels of all the wheels whose speeds are detected are lowered. 

In the light of the difficulty in detecting a change in the tire pressure of a vehicle wheel, J P-A-62- 149502 proposes 
a tire pressure detecting apparatus which operates depending on a fact that the tire pressure changes with an amount 
or rate of change of the rotating speed of the corresponding wheel when an external force acting on the wheel varies 
for some reason, for example, when the wheel runs over a projection of the road surface. Since the maximal value of 

35 the amount of change of the rotating speed of the wheel decreases with a decrease in the tire pressure, the tire pressure 
can be estimated on the basis of the detected maximal value of the wheel speed change amount. The tire pressure 
detecting apparatus includes signal generating means for generating a signal when the detected amount of change of 
the wheel speed is higher than a predetermined threshold value, and tire pressure estimating means responsive to that 
signal, for estimating the tire pressure on the basis of the maximal value of the wheel speed change amount. 

40 However, the condition of the road surface that permits the above-indicated apparatus to detect the tire pressure is 
limited. Described in detail, the tire pressure can be estimated from the maximal value of the detected wheel speed 
change amount, only where the relevant wheel runs over a single projection on the road surface. When the wheel passes 
successive projections and recesses on a bumpy or washboard road surface, changes of the rotating speed of the wheel 
due to those projections and recesses overlap each other, whereby the maximum value of the detected wheel speed 

45 change amount is not commensurate with the tire pressure, making it impossible to effect accurate detection of the tire 
pressure. 

While the difficulty of detecting the tire pressure of the vehicle wheel has been described above by way of example, 
it has been also difficult to detect other parameters indicative of the motion of the vehicle wheel, such as the angular 
velocity, angular acceleration, vertical or normal speed, and vertical or normal acceleration of the wheel, and the changing 

so states of the other characteristics of the vehicle wheel and/or tire, such as the eccentricity of the wheel, replacement of 
the metal wheel member per se, wear of the tire, and attachment of a foreign matter to the tire. 

In view of the above situation, the present inventors proposed a wheel information estimating apparatus which 
includes (a) a detecting device for detecting a wheel motion state or parameter indicative of the motion of a tired wheel 
of a motor vehicle; (b) a disturbance observer for estimating a disturbance acting on the wheel, on the basis of at least 

55 the detected wheel motion state, and a basic value indicative of a specific kind of wheel information about the tired 
wheel; (c) means for estimating an amount of change from the basic value to a current value indicative of a current state 
of the wheel information of the tired wheel, on the basis of the disturbance detected by the disturbance observer, and 
(d) means for estimating the wheel information which is a sum of the estimated amount of change and the basic value 
indicated above. 



2 



BNSDOCID <EP 0699546A1_I_> 



EP 0 699 546 A1 

The disturbance which acts on the tired wheel is not limited to that given to the wheel from the road surface due to 
its irregularity. That is, the tired wheel may also suffer from other disturbance components such that those due to a 
change or variation in the moment of inertia or spring constant of the wheel. The moment of inertia of the tired wheel 
varies due to wear of the tire, attachment of some foreign matter to the tire, replacement of the tire or metal wheel 
5 member with a new one, and the like. The spring constant of the tired wheel varies due to a change of the air pressure 
of the tire and the like. 

As a result of a further study on the wheel information estimating apparatus as described above, the present inventors 
found the following facts. That is, the disturbance observer is constructed using a dynamic model of a tired wheel. If the 
dynamic model is required to represent a lot of vibrational modes relating to the wheel, the construction of the disturbance 

10 observer becomes complicated, making it difficult to improve its estimating accuracy. On the other hand, the wheel 
motion state has a frequency characteristic which has a certain relationship with the estimating accuracy of the distur- 
bance observer for estimating the specific kind of information about the tired wheel. The present inventors found that if 
the frequency characteristic of the wheel motion state is optimally controlled, the estimating accuracy of the disturbance 
observer can be easily improved, without using such a complicated dynamic model of the wheel. 

75 Another type of wheel information estimating apparatus is proposed in JP-A-5-133831 . This apparatus includes a 
detecting device adapted to detect a rotating speed of a tired wheel as a wheel motion state, and estimating means for 
estimating an air pressure of a tire as a specific kind of wheel information, on the basis of a frequency characteristic of 
the wheel motion state, which is received from the detecting device. More specifically, the tire air pressure is estimated 
on the basis of the frequency of a selected one of a plurality of frequency components of the wheel motion state, which 

20 selected one having substantially the greatest strength within a predetermined frequency range. The thus constructed 
apparatus is referred to as a resonance-frequency type wheel information estimating apparatus. 

The present inventors then found that the above-described technique of optimally controlling the frequency charac- 
teristic may be effectively applied to the wheel information estimating apparatus of the type as described just above. 
According to the findings of the inventors, a filter or other device may be interposed between the detecting device and 

25 the estimating means, so as to modify the frequency characteristic of the wheel motion state, by particularly emphasizing 
the selected one of the plural frequency components of the parameter, which is within the predetermined frequency 
range. The parameter having the thus modified frequency characteristic is applied to the estimating means, with a 
reduced burden on the estimating means, which leads to improved estimating accuracy. 

That is, with respect to the wheel information estimating apparatus which generally includes (a) a detecting device^ 

30 for detecting a wheel motion state, and (b) an estimating device for estimating a specific wheel information on the basis 
of the detected wheel motion state, the present inventors found it possible to employ the above technique of optimally 
controlling the frequency characteristic, so as to improve the estimating accuracy of the apparatus, no matter whether 
the apparatus is of the former type having the disturbance observer, or of the latter type, namely, of the resonance- 
frequency type. 

35 In the former type of wheel information estimating apparatus having the disturbance observer, the basic value indic- 
ative of the specific wheel information is fixed or predetermined, and the disturbance observer estimates the disturbance 
on the basis of the fixed basic value. A further study by the present inventors on this type of apparatus revealed that the 
accuracy of estimating the wheel-information can be improved by changing or compensating the wheel information basic 
value so that the amount of change from the basic value to the current value of the wheel information is zeroed, whereby 

40 the basic value used by the disturbance observer accurately represents an actual or current state of the wheel informa- 
tion. This method is also effective to improve the estimating accuracy of the wheel information estimating apparatus of 
the former type. 

SUMMARY OF THE INVENTION 

45 

It is therefore an object of the present invention to provide a wheel information estimating apparatus which is able 
to estimate a specific kind of information relating to a tired wheel of a motor vehicle, with significantly improved accuracy. 

The above object may be accomplished according to the principle of the present invention, which provides a wheel 
information estimating apparatus, comprising: a detecting device for detecting a wheel motion state relating to a move- 
so ment of a tired wheel of a motor vehicle; estimating means for estimating wheel information relating to the tired wheel, 
on the basis of the wheel motion state detected by the detecting device; and modifying means for modifying an element 
relating to estimation of the wheel information, so as to improve an accuracy of estimating the wheel information, the 
element being at least one of: a movement of a center of rotation of the tired wheel as viewed in a plane in which the 
wheel is rotated; a signal applied from the detecting device to the estimating means; and an internal parameter used by 
55 the estimating means for estimating the wheel information. 

The wheel motion state as indicated above may be selected from among an angular velocity, angular acceleration, 
vertical velocity, vertical acceleration, longitudinal velocity, and longitudinal acceleration. With respect to a tired wheel 
which has a rim and a belt that are connected to each other through a torsion spring, for example, the wheel information 
as indicated above may be selected from among: an air pressure of a tire of the tired wheel; a radius of the tire; the 
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moment of inertia of the tired wheel, more specifically, the moment of inertia of the rim and that of the belt; a spring 
constant of the torsion spring; a damping coefficient of the tire; a degree of contact of the tired wheel with a road surface, 
and a cornering power of the tire. 

In one advantageous form of the wheel information estimating device as described above, the modifying means 
includes frequency characteristic modifying means for providing an optimum frequency characteristic of the wheel motion 
state to be supplied to the estimating means, by increasing a ratio of a strength of a selected one of a plurality of frequency 
components of the wheel motion state, the selected one being within a selected frequency range, to a strength of the 
other frequency components outside the selected frequency range. In the above form of the invention, the frequency 
characteristic modifying means selects one of the plural frequency components of the wheel motion state, which is within 
the selected frequency range, increases the ratio of the strength of the selected frequency component to that of the 
other frequency components outside the selected frequency range, and applies the thus modified frequency character- 
istic to the estimating means. Accordingly, the accuracy of estimating the wheel information can be easily improved by 
setting the above selected frequency range to a frequency range in which the estimating means is able to estimate the 
wheel information with high accuracy. Since the estimating means receives the optimum frequency characteristic of the 
wheel motion state, the present apparatus is able to utilize a simplified dynamic model of the tired wheel for estimation 
of the wheel information, thus assuring sufficiently high estimating accuracy. The above-indicated ratio of the strength 
of the selected frequency component of the wheel motion state to that of the other components may be increased either 
by increasing the strength of the selected frequency component, or reducing the strength of the other frequency com- 
ponents. 

In the above form of the present invention, the frequency characteristic modifying means may include means for 
changing the selected frequency range, on the basis of at least one of an estimated value of the wheel information 
obtained by the estimating means, and a related value of the estimated value. 

The frequency characteristic representing the wheel motion state is not always constant, but varies with a vehicle 
speed, for example. If the selected frequency range is fixed regardless of such variation of the frequency characteristic, 
therefore, an optimum frequency characteristic suitable for estimation of the wheel information cannot be obtained. In 
view of this, the selected frequency range is changed by the above-indicated means, on the basis of the estimated value 
obtained by the estimating means, which reflects a change in the frequency characteristic of the wheel motion state. 
Thus, the selected frequency range varies with a change in the actual frequency characteristic of the wheel motion state. 

In the above form of the present invention, the frequency characteristic modifying means may also include means 
for changing the selected frequency range, so that an amount of variation of an estimated value of the wheel information 
obtained by the estimating means is not greater than a predetermined value. That is, the selected frequency range is 
changed in a feed-back fashion so as to reduce the amount of change of the estimated value. Thus, the selected fre- 
quency range is always suitably determined to be suited with the actual frequency characteristic of the wheel motion state. 

The frequency characteristic of the wheel motion state has a certain relationship with a locus of the center of a tired 
wheel as seen in the transverse direction of a motor vehicle, which is drawn upon relative movements of the tired wheel 
and a vehicle body through a suspension. Where a wheel speed that is a rotating speed of the wheel is chosen as the 
wheel motion state, and the above locus of the center of the wheel includes a front-to-rear or longitudinal component 
extending in a longitudinal or running direction of the vehicle as well as a vertical component, the wheel speed varies 
due to vibrations in the longitudinal direction of the vehicle, which are applied to the outer peripheral surface of the wheel 
due to wheel vibrations caused by projections and recesses of the road surface. As a result, the influence of the wheel 
vibrations appear in a particular frequency range of the wheel speed. Accordingly, if the geometry and the like of the 
suspension is designed such that the center of a circle (e.g., instantaneous circle) which approximates the locus of the 
center of the wheel is located as close as possible to the center of the wheel, so that the locus of the center of the wheel 
consists as much as possible of the front-to-rear or longitudinal component, the influence of the above wheel vibrations 
on the wheel speed is relatively emphasized in the particular frequency range. Conversely, if the geometry of the sus- 
pension is designed such that the center of the above- indicated approximate circle is located as far as possible from 
the center of the wheel, so that the locus of the center of the wheel consists as much as possible of the vertical component, 
the influence of the wheel vibrations on the wheel speed is relatively minimized in the particular frequency range. 

The above-indicated frequency characteristic modifying means may generally consist of a digital filter or analog 
filter. In this case, such a filter is disposed between the detecting device and the estimating means. Since there is a 
certain relationship between the locus of the center of the wheel and the frequency characteristic of the wheel motion 
state, as described above, an optimum frequency characteristic may be achieved by appropriately designing or otherwise 
adjusting the suspension itself, without using the filter as described above. 

In the above form of the wheel information estimating apparatus of the present invention, therefore, the frequency 
characteristic modifying means may include a geometry of a suspension which connects the tired wheel and a vehicle 
body such that the wheel and the vehicle body are displaceable relative to each other, the geometry being designed 
such that a center of the tired wheel provides an optimum locus as seen in a transverse direction of the motor vehicle, 
so as to improve accuracy of estimating the wheel information. The optimum locus of the center of the tired wheel, which 
is drawn upon a relative movement between the wheel and the vehicle body, results in an optimum frequency charac- 
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teristic of the wheel motion state. This eliminates a need to provide such a filter as described above, or other device or 
program for providing the optimum frequency characteristic. 

The wheel information includes some kinds of information which can be estimated with effectively improved accuracy, 
by relatively emphasizing the frequency component of the wheel speed (wheel motion state) which is most influenced 
5 by the wheel vibrations. When the estimating means estimates these kinds of information, the geometry of the suspen- 
sion, for example, can be designed so that the influence of the wheel vibrations on the wheel speed appears particularly 
strongly in the frequency range, that is, so that the locus of the center of the wheel consists as much as possible of the 
front-to-rear or longitudinal component. This arrangement is one example of the above-indicated suspension adjusting 
means. 

w The wheel information includes other kinds of information which can be estimated with effectively improved accuracy, 
by relatively weakening the frequency component of the wheel speed (wheel motion state) which is most influenced by 
the wheel vibrations. When the estimating means estimates these kinds of information, the geometry of the suspension, 
for example, can be designed so that the influence of the wheel vibrations on the wheel speed hardly appears in the 
frequency range, that is, so that the locus of the center of the wheel consists as much as possible of the vertical com- 

15 ponent, and includes almost no front-to-rear component. This arrangement is another example of the above-indicated 
suspension adjusting means. 

The latter example of the suspension adjusting means may provide a suspension which includes an upper arm and 
a lower arm, which cooperate to support a tired wheel relative to a vehicle body, such that the wheel is displaceable 
relative to the vehicle body. The geometry of this suspension may be designed such that the upper and lower arms 

20 extend substantially in parallel to each other as seen in the transverse direction of the vehicle, so that the locus of the 
center of the wheel does not substantially include a longitudinal or front-to-rear component. 

In the above form of the wheel information estimating device of the present invention, the estimating means may 
include (a) a disturbance observer for estimating a disturbance acting on the tired wheel, on the basis of a basic value 
of the wheel information, and the wheel motion state received from the detecting device, and (b) means for estimating 

25 an amount of change from the basic value to a value indicative of an actual state of the wheel information, on the basis 
of the disturbance estimated by the disturbance observer. Since the disturbance observer receives the wheel motion 
state whose frequency characteristic is appropriately controlled by the frequency characteristic modifying means, the 
disturbance observer is able to estimate the disturbance with high accuracy, whereby the estimating means can always 
estimate the wheel information with high accuracy. 

30 In the above form of the present invention, the estimating means may estimate the wheel information, on the basis 
of a frequency of a selected one of a plurality of frequency components of the wheel motion state received from the 
detecting device, the selected one of the frequency components having substantially the greatest strength within a 
selected frequency range. Since the estimating means receives the wheel motion state whose frequency characteristic 
is appropriately controlled by the frequency characteristic modifying means, the estimating means is able to obtain an 

35 appropriate frequency needed for estimation with constantly high accuracy, which leads to improved accuracy of esti- 
mating the wheel information. 

In the above form of the present invention, the detecting device comprises a wheel speed detector for detecting, as 
the wheel motion state, a wheel speed which is a rotating speed of the tired wheel, and supplying the wheel speed to 
the estimating means, the estimating means estimating a degree of contact of the tired wheel with a road surface as 

40 the wheel information, on the basis of a strength of a selected one of a plurality of frequency components of the wheel 
speed received from the wheel speed detector, the selected one of the frequency components being within a selected 
frequency range. It is noted that the strength of the frequency component of the wheel speed which is within the selected 
frequency range increases with reduction in the degree of contact of the tired wheel with the road surface, that is, the 
road holding characteristic of the wheel. Using this relationship, the degree of contact of the wheel with the road surface 

45 is estimated on the basis of the detected wheel speed. Since the estimating means receives the wheel speed whose 
frequency characteristic is appropriately controlled by the frequency characteristic modifying means, the estimating 
means is able to obtain the strength of the frequency component needed for the estimation, which leads to improved 
accuracy of estimating the degree of contact between the wheel and the road surface. 

In the above-described arrangements of the above form of the invention as described above, only one estimated 

so value obtained by the estimating means is used as a final value to be generated by the same means. In this instance, 
even if the frequency characteristic of the wheel motion state is appropriately controlled as described above, the esti- 
mated value may vary from one cycle to another cycle of the routine for estimating the wheel information, resulting in 
reduced estimating accuracy. In view of this situation, the present wheel information estimating apparatus may be 
adapted such that the estimating means obtains a plurality of provisional estimated values of at least one of the wheel 

55 information and a parameter relating to the wheel information, so as to obtain a current final estimated value of the wheel 
information, on the basis of the plurality of provisional estimated values. A certain number of the provisional estimated 
values are used for at least one of re-estimation of the currentf inal estimated value and estimation of a next final estimated 
value of the wheel information. The number of the provisional estimated values is variable depending upon whether an 
amount of change from a previously estimated value of the wheel information to a currently estimated value thereof is 
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larger than a predetermined value. Since only one final value is obtained from a plurality of provisional estimated values 
of the wheel information, the obtained final value is less likely to be influenced by the variation in the individual provisional 
es mated values. Further, the number of the provisional values used for obtaining the one final value is properly changed 
as needed, depending upon an error in estimating the provisional values, so as to assure sufficiently high estimating 
ac jracy of the apparatus. 

In one arrangement of the wheel information estimating apparatus as described just above, the estimating means 
may be adapted to determine the number of provisional values used for obtaining the final value, depending upon the 
amount of change of the currently estimated value. In another arrangement, the estimating means may be adapted to 
increase the number of provisional values in increment, until the amount of change of the currently estimated value 
becomes equal to or smaller than the predetermined value, or increase the number of provisional values depending 
upon a difference between the amount of change and the threshold value. 

In another form of the wheel information estimating apparatus of the present invention, the modifying means modifies 
the internal parameter used by the estimating means, and the estimating means includes (a) a disturbance observer for 
estimating a disturbance acting on the tired wheel, on the basis of at least the wheel motion state detected by the 
detecting device, and a basic value of the wheel information, the disturbance observer using the basic value as the 
internal parameter, (b) change amount estimating means for estimating an amount of change from the basic value to a 
value indicative of an actual state of the wheel information, on the basis of the disturbance estimated by the disturbance 
observer, (c) basic value compensating means as the modifying means for compensating the basic value of the wheel 
information to provide a compensated wheel information basic value, on the basis of the amount of change estimated 
by the change amount estimating means, so that the estimated amount of change is not larger than a predetermined 
value; and (d) wheel information estimating means for estimating the compensated wheel information basic value as 
the wheel information. 

Where the basic value of the wheel information is determined as a fixed value, the accuracy of estimating an amount 
of change of the wheel information is reduced with an increase in the estimated amount of change from the basic value 
to the currently estimated value indicative of a current or actual state of the wheel information. This results in reduced 
accuracy of estimating the wheel information. If the wheel information basic value is changed so as to be close to the 
current value, thereby to reduce the estimated amount of change to a minimum, the wheel information can be estimated 
with sufficiently high accuracy. In the above form of the invention, therefore, the basic value compensating means is 
adapted to compensate the wheel information basic value to be used by the disturbance observer, on the basis of the 
estimated amount of change of the wheel information, so that the estimated change amount becomes equal to or smaller 
than the predetermined value, thereby to provide a compensated wheel characteristic basic value. The estimating means 
then estimates the wheel information on the basis of the thus compensated basic value. In this arrangement, the dis- 
turbance observer is able to estimate a disturbance with improved accuracy, thereby assuring improved accuracy of 
estimating the wheel information. 

The disturbance observer estimates a disturbance as one of variables of state. The disturbance observer may 
assume a model of the tired wheel to be an integral rotary body on which a disturbance acts, or a model having a rim 
and a belt which are coupled with each other through a torsion spring, such that the rim and belt are rotatable relative 
to each other. The disturbance observer may also assume a model of the wheel having a rim and a belt which are 
coupled with each other, through a torsion spring and a damper which are connected in parallel with each other, such 
that the rim and belt are rotatable relative to each other. The disturbance observer may also assume a disturbance to 
act on a tired wheel in which a wheel (consisting of a disc and a rim) and a tire (consisting of a belt and side walls) are 
coupled with each other by a spring, such that the wheel and the tire are displaceable relative to each other in the vertical 
direction of the vehicle. In another model which may be assumed by the disturbance observer, a disturbance acts on a 
tired wheel in which a wheel and a tire are coupled with each other, by a torsion spring and a damper which are connected 
in parallel with each other, such that the wheel and the tire are displaceable relative to each other in the vertical direction 
of the vehicle. 

When the disturbance observer assumes a model of the tired wheel having a rim and a belt which are connected 
to each other by a torsion spring, such that the rim and the belt are rotatable relative to each other, the disturbance 
observer may estimate an angular velocity co B of the belt, a torsion angle 6 RB between the rim and the belt, and a 
disturbance w 2 . on the basis of a detected value of an angular velocity o>r of the rim as the wheel motion state, and 
basic values of a moment of inertia J R of the rim, a moment of inertia J B of the belt, and a spring constant K of the torsion 
spring, as the wheel information basic value. In this case, the change amount estimating means estimates an amount 
of change AK from the basic value of the spring constant K to a value indicative of an actual spring constant, on the 
basis of estimated values of the angular velocity co B of the belt, the torsion angle 8 RB , and the disturbance w 2 . and the 
basic value compensating means compensates the basic value of the spring constant K, so that the estimated amount 
of change AK is not larger than a predetermined value. Then, the wheel information estimating means estimates the 
basic value of the spring constant K as the wheel information, when the amount of change AK is not larger than the 
predetermined value. 
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When the disturbance observer assumes an equivalent model of the tired wheel which is equivalent to a basic model 
in which a rim is coupled with a belt through a torsion spring and a damper which are connected in parallel with each 
other, such that the rim and the belt are rotatable relative to each other, the equivalent model being designed such that 
a rim is coupled with a belt through a torsion spring and a damper which are connected in parallel with each other, such 

5 that the rim and the belt are displaceable relative to each other in one linear direction, the disturbance observer may 
estimates an equivalent linear displacement velocity x B of the belt, an equivalent relative linear displacement x RB 
between the rim and the belt, and a disturbance w 2 . on the basis of a detected value of an equivalent linear displacement 
velocity x R * of the rim, as the wheel motion state, and basic values of an equivalent inertial mass m R of the rim, an 
equivalent inertial mass m B of the belt, a spring constant K w and a damping coefficient D w , as the wheel information 

jo basic value. In this case, the change amount estimating means estimates an amount of change AK W from the basic 
value of the spring constant K w to a value indicative of an actual spring constant, and an amount of change ADw from 
the basic value of the damping coefficient D w to a value indicative of an actual damping coefficient, on the basis of 
estimated values of the equivalent linear displacement velocity x B of the belt, the equivalent relative linear displacement 
x RB and the disturbance w 2 . and the basic value compensating means compensates the basic values of the spring 

is constant K w and the damping coefficient D^ so that the estimated amounts of change AK W and AD W are not larger 
than respective predetermined values. Then, the wheel information estimating means estimates the basic values of the 
spring constant K w and the damping coefficient D w as the actual spring constant and damping coefficient, when the 
estimated amounts of change AK W and AD W are not larger than the respective predetermined values. 

20 BRIEF DESCRIPTION OF THE DRAWINGS 

The above and optional objections, features and advantages of the present invention will be better understood by 
reading the following detailed description of presently preferred embodiments of the invention, when considered in con- 
nection with the accompanying drawings, in which: 

25 

Fig. 1 is a schematic block diagram illustrating one embodiment of a wheel information estimating apparatus of the 
present invention adapted to estimate certain wheel information relating to a certain condition of tired wheels of an 
automotive vehicle; 

Fig . 2 is a fragmentary cross sectional view of a tired vehicle wheel whose disturbances are detected by the apparatus 
30 of Fig. 1 ; 

Fig. 3 is a view showing a dynamic model of the tired wheel of Fig. 2; 

Fig. 4 is a block diagram illustrating the functional elements of the apparatus of Fig. 1 ; 

Fig. 5 is a flow chart illustrating a routine for obtaining a rotating speed variation of the tired wheel according to a 
control program stored in a read-only memory of a computer used by the apparatus; 
35 Fig. 6 is a flow chart illustrating a wheel speed calculating/compensating routine according to another control pro- 
gram also stored in the read-only memory; 

Fig. 7 is a graph illustrating a relationship between the number N used in step S8 of the routine of Fig. 5 and vehicle 
speed V; 

Fig. 8 is a graph explaining approximation of the disturbances as detected by the apparatus; 
40 Fig. 9 is a block diagram illustrating a disturbance observer incorporated in the apparatus; 

Fig. 10 is a flow chart illustrating a tire pressure alarming routine according to a control program stored in a read- 
only memory of another computer used by the apparatus of Fig. 1 ; 

Fig. 11 is a flow chart illustrating a routine for calculating correlation functions for obtaining a spring constant variation 
of the tired wheel, according to another control program stored in the read-only memory; 
45 Fig. 12 is a graph illustrating a frequency characteristic of a wheel speed signal received by a pre-processing filter 

of the apparatus of Fig. 1 ; 

Fig. 13 is a graph illustrating a frequency characteristic of a wheel speed signal generated as an output from the 
pre-processing filter; 

Fig. 14 is a graph illustrating a filtering characteristic of the pre-processing filter; 
so Fig. 1 5 is a graph illustrating a relationship between a variation 5 as calculated in step S40 of the routine of Fig. 10, 

and the reference value "M"; 

Fig. 1 6 is a blockdiagram illustrating the functional elements of a wheel information estimating apparatus constructed 
according to a second embodiment of the invention; 

Fig. 1 7 is a block diagram illustrating the construction of the wheel information estimating apparatus of Fig. 1 6; 
55 Fig. 18 is a flow chart illustrating a tire pressure alarming routine according to a control program stored in a read- 
only memory of a computer used by the apparatus of Fig. 16; 

Fig. 19 is a graph showing a relationship between a variation 5 as calculated in step S60 of the routine of Fig. 18, 
and the reference value N; 
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Fig. 20 is a flow chart illustrating another tire pressure alarming routine according to a control program stored in a 
read-only memory of a computer, which is used by a third embodiment of a wheel information estimating apparatus 
of the present invention; 

Fig. 21 is a graph showing a relationship between an amount of change Afo of the resonance frequency and the 
5 number "N" of final values to be obtained in the routine of Fig. 20; 

Fig. 22 is a block diagram illustrating the functional elements of a wheel information estimating apparatus constructed 
according to a fourth embodiment of the invention; 

Fig. 23 is a graph explaining the function of a pre-processing filter of the apparatus of Fig. 22; 
Fig. 24 is a flow chart illustrating a tire pressure alarming routine according to a control program stored in a read- 
io only memory of a computer used by the apparatus of Fig. 22; 

Fig. 25 is a flow chart illustrating a tire pressure alarming routine according to another control program stored in the 
read-only memory; 

Fig. 26 is a flow chart illustrating a tire pressure alarming routine according to a further control program stored in 
the read-only memory; 

15 Fig. 27 is a perspective view showing a rear suspension system of an automotive vehicle; 

Fig. 28A is a view showing an instantaneous center O located far from the center of a vehicle wheel, and a locus 
of the center of the wheel, Fig. 28B showing an instantaneous center O located in a normal position, and a locus 
of the center of the wheel, Fig. 28C showing an instantaneous center O located closer to the center of the wheel, 
and a locus of the center of the wheel; 

20 Fig. 29 is a perspective view showing an upper arm and lower arms of a suspension which provides a substantially 

straight locus of the center of the wheel; 

Fig. 30 is a perspective view showing an upper arm and lower arms of a suspension which provides a normally 
curved locus of the center of the wheel; 

Fig. 31 is a view showing three kinds of locuses of the center of the wheel, which correspond to those of Figs. 28A, 
25 28B and 28C; 

Fig. 32 is a graph illustrating a relationship between each locus of the center of the wheel as shown in Fig. 31 , and 
a frequency characteristic of a wheel speed signal for each wheel; 

Fig. 33 is a block diagram illustrating the functional elements of a wheel information estimating apparatus constructed 
according to a fifth embodiment of the invention; 
30 Rg. 34 is a flow chart illustrating a routine for calculating correlation functions for obtaining variation of a moment 
of inertia of a belt of a tired wheel, according to a control program stored in a read-only memory of the apparatus 
of Fig. 33; 

Fig. 35 is a flow chart illustrating a routine for calculating correlation functions for obtaining variation of a moment 
of inertia of a rim of the wheel, according to another control program also stored in the read-only memory; 
35 Fig. 36 is a block diagram illustrating the functional elements of a wheel information estimating apparatus constructed 
according to a sixth embodiment of the invention; 

Fig. 37 is a view showing a tire-suspension model used by a disturbance observer of the apparatus of Fig. 36; 
Fig. 38 is a graph showing frequency characteristics of wheel speed signals received by a pre-processing filter of 
Fig. 36 before and after a change in the moment of inertia, respectively; 
40 Fig. 39 is a graph illustrating frequency characteristics of the wheel speed signals generated as outputs from the 
pre-processing filter; 

Fig. 40 is a blockdiagram illustrating the functional elements of a wheel information estimating apparatus constructed 
according to a seventh embodiment of the invention; 

Fig. 41 is a graph illustrating a frequency characteristic of a wheel speed signal obtained when a tired wheel is not 
45 in a good contact with a road surface; 

Fig. 42 is a graph illustrating a frequency characteristic of a wheel speed signal obtained when the tired wheel is in 
a good contact with the road surface; 

Fig. 43 is a flow chart illustrating a routine for determining a degree of contact of the tired wheel with the road surface, 
according to a control program stored in a read-only memory of a computer used by the apparatus of Fig. 40; 
so Fig. 44 is a schematic blockdiagram illustrating an eighth embodiment of a wheel information estimating apparatus 

of the present invention; 

Fig. 45 is a blockdiagram illustrating the functional elements of the wheel information estimating apparatus of Fig. 44; 
Fig. 46 is a flow chart illustrating a tire pressure detecting routine according to a control program stored in a read- 
only memory of a computer used by the apparatus of Fig. 44; 
55 Fig. 47 is a view showing a dynamic model of a disturbance observer incorporated in a ninth embodiment of a wheel 
information estimating apparatus of the invention; 
Fig. 48 is a simplified view of the dynamic model of Fig. 47; 

Fig. 49 is a blockdiagram showing the functional elements of the apparatus according to the embodiment of Fig. 47; 
Fig. 50 is a block diagram depicting a computer used in the apparatus of Fig. 49; 
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Fig. 51 is a flow chart illustrating a routine for detecting air pressure of the tired wheel according to a control program 
stored in the read-only memory of the computer of Fig. 50; and 

Fig. 52 is a flow chart illustrating a modified tire pressure detecting routine according to a control program stored in 
the read-only memory. 

5 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Referring first to Figs. 1 and 2, the wheel information estimating apparatus constructed according to the first embod- 
iment of this invention includes a magnetic pickup (variable-reluctance pickup) 1 2 disposed adjacent to a rotor 10 which 

w rotates with a tired wheel indicated generally at 14 in Fig. 2. The rotor 10 has a multiplicity of teeth 16 formed along its 
outer periphery such that the teeth 16 are equally spaced apart from each other in the circumferential direction, that is, 
in the direction of rotation of the rotor 10. 

The magnetic pickup 12 is constructed to generate a voltage signal whose amplitude periodically changes as the 
teeth 16 pass the sensing head of the pickup 12. The waveform of the voltage signal generated by the pickup 12 is 

75 shaped into a rectangular pulse form by a waveform shaper 1 8, and the shaped pulse signal is applied to an input/output 
port 22 of a computer 20. While only one set of the rotor 10, pickup 12 and waveform shaper 18 is shown in Fig. 1, for 
illustrative purpose only, the motor vehicle has four tired wheels 14, and the rotor 10, pickup 12 and waveform shaper 
1 8 are provided for each of the four tired wheels 14. Each of the four pickups 12 is connected to the computer 20 through 
the appropriate waveform shaper 1 8. 

20 Each tired wheel 14 consists of a metal wheel member 24, and a tire 26 attached to a rim of the wheel member 24. 
The tired wheel 14 is dynamically simulated as a system or model as illustrated in Fig. 3, in which a rim 28 and a belt 
30 are connected to each other by a torsion spring 32 such that the rim 28 and the belt 30 are rotatable relative to each 
other. The belt 30 is considered to be the outer surface of the tire 26 at which the tired wheel 14 contacts the road 
surface. Since the rotor 10 is disposed so as to rotate with the metal wheel member 24, the magnetic pickup 12 is 

25 adapted to detect an angular velocity of the rim 28, in a strict sense. 

The computer 20 incorporates a CPU 40 as a central processing unit, a ROM 42 in the form of a read-only memory 
as a first memory device, and a RAM 44 in the form of a random-access memory as a second memory device, as shown 
in Fig. 1. The ROM 42 stores control programs for executing the routines illustrated in the flow charts of Figs. 5 and 6, 
and cooperates with the CPU 40 and RAM 44 to constitute a rim speed calculator/compensator 45 indicated in the 

30 schematic block diagram of Fig. 4. 

The computer 20 is connected to another computer 47 as shown in Fig. 1 , which incorporates a CPU 48 as a central 
processing unit, a ROM 49 in the form of a read-only memory as the first memory device, a RAM 50 in the form of a 
random-access memory as the second memory device, and an input/output port 51 as an input/output device. The ROM 
49 stores various control programs for executing the routines, such as a tire pressure alarming routine as illustrated in 

35 the flow chart of Fig. 10, and cooperates with the CPU 48 and RAM 50 to constitute a disturbance observer 52, a 
parameter calculator 53 including a correlation calculator 56 and a normalizer 58, and a determining processor 62, as 
indicated in the block diagram of Fig. 4. 

The computers 20, 47 are connected to each other through a pre-processing filter 64. More specifically described 
referring to Fig. 4, the rim speed calculator/compensator 45 is connected to the disturbance observer 52 through the 

40 pre-processing filter 64, such that an output signal indicative of a rotating speed of each tired wheel 14 is transmitted 
from the rim speed calculator/compensator 45 to the disturbance observer 52, through the pre-processing filter 64. While 
the pre-processing filter 64 principally consists of a digital signal processor (DSP) in the present embodiment, the filter 
64 may principally consist of a discrete circuit. The operation of the pre-processing filter 64 will be described in detail later. 
The input/output port 51 of the computer 47 is connected to a display device 66 which provides information indicative 

45 of a result of judgement or determination made by the determining processor 62, which would aid an operator of the 
vehicle in recognizing the condition of each tired wheel 14. While the display device 66 is a liquid crystal display in the 
present embodiment, the liquid crystal display device 66 may be replaced by any other suitable indicator or display such 
as an indicator light device operated in a predetermined manner, for instance, in a flickering manner, and a voice indicator 
which speaks to the vehicle operator. 

so The input/output port 51 is also connected to a torque detector 68 adapted to detect a driving torque or a braking 
torque of the wheel member 24 (rim 28). For example, the torque detector 68 may include a strain gage fixed to the axle 
of the wheel member 24. 

The rim speed calculator/compensator 45 operates to calculate the rotating speed of each of the four tired wheels 
14. on the basis of the signals received from the corresponding magnetic pickups 12 through the respective waveform 
55 shapers 1 8. The calculator/compensator 45 is further adapted to adjust the calculated speed of each tired wheel 1 4, on 
the basis of the vehicle running speed (speed of a body 72 of the vehicle as indicated in Fig. 4) detected by a suitable 
vehicle speed sensor 70 as indicated in Fig. 1. The vehicle speed sensor 70 may be a Doppler-effect ground-speed 
sensor. The adjustment of the speed of the tired wheel 14 as calculated by the magnetic pickup 12 is effected to com- 
pensate the calculated speed for a periodic or cyclic variation of the wheel speed which would arise from inherent 
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deviations of the tired wheel 14 and rotor 10 from the nominal values due to errors during manufacture and assembling 
thereof. 

The rotating speed of the tired wheel 14 or wheel speed is obtained as the peripheral speed of the tire 26. To this 
end, the radius of the tire 26 (i.e.. the distance from the road surface to the axis or center of the wheel 14) should be 

5 known. In this respect, the radius of the tire 26 varies with the amount of deformation of the tire due to a load acting 
thereon, which in turn varies with the air pressure in the tire 26. As long as the air pressure of the tire 26 remains at the 
nominal or desired level, the nominal radius of the tire 26 is used to calculate the peripheral speed of the tired wheel 1 4 
as the rotating speed. If a change of the air pressure larger than a predetermined limit is found in an appropriate routine 
as described later, the radius of the tire 26 used to calculate the peripheral speed is determined according to a prede- 

10 termined relationship (stored in the ROM 42) between the tire radius and the amount of change of the air pressure of 
the tire. 

The rim speed calculator/compensator 45 is adapted to execute a routine for obtaining the periodic wheel speed 
variation of Fig. 5 and a routine for calculating and compensating the detected wheel speed of Fig. 6. 

The wheel speed variation obtaining routine of Fig. 5 is executed at least once after the rotor 10 and the wheel 14 

15 have been assembled on the vehicle. An operation to obtain the wheel speed variation according to this routine may be 
performed by a manufacturer of the vehicle, a service engineer of a repair shop, or a user of the vehicle. Where the 
wheel speed variation is obtained by the user, the routine is executed each time a predetermined condition is satisfied, 
for example, each time the cumulative running distance or time of the vehicle reaches a predetermined value. In any 
case, it is desirable to execute the routine while the vehicle is running in a stable mode at a substantially constant speed 

20 within a predetermined range, without acceleration or deceleration (braking). 

The wheel speed variation obtaining routine of Fig. 5 will be first described. The routine is started with step S1 in 
which values "n", V, "V" and "v n " are initialized. Step S1 is followed by step S2 to calculate the vehicle speed V on the 
basis of the output of the vehicle speed sensor 70. Then, the control flow goes to step S3 to calculate the speed v n of 
the tired wheel 14 (hereinafter referred to as "wheel speed v n "). In a strict sense, the wheel speed v n calculated is the 

25 peripheral speed of the belt 30 (peripheral speed of the tire 26 at its outer circumference contacting the road surface), 
assuming that an angular velocity of the belt 30 is the same as an angular velocity o> R of the rim 28. For easier 
understanding, it is assumed that the wheel speed v n is calculated when each one of the teeth 16 passes the sensing 
head of the magnetic pickup 12. More specifically described, the wheel speed v n corresponding to each of the teeth 16 
of the rotor 10 is calculated, on the basis of a time interval between the rises or falls of the adjacent two rectangular 

30 pulses obtained from the waveform shaper 18, or a time interval between midpoints between the rise and fall of the 
adjacent pulses. 

However, the relation between the time required for each tooth 16 to pass the pickup 12 and the time required to 
calculate the wheel speed v n should be taken into account. In this respect, the routine of Fig. 5 may be modified such 
that the wheel speed v n is calculated each time two or more teeth 16 pass the pickup 12. 

35 Step S3 is followed by step S4 to calculate a difference (V - v n ) of the calculated wheel speed v n from the vehicle 
speed V determined in step S2. The calculated speed difference (V - v n ) is considered a variation of the wheel speed 
v n (hereinafter referred to as "wheel speed variation (V - v n )"). The thus obtained wheel speed variation (V - v n ) is stored 
at the appropriate one of successive addresses of a "wheel speed variation" memory of the RAM 44 whose number is 
equal to T (number of the wheel speed values v n obtained during one revolution of the rotor 10 or wheel 14). In the 

40 present embodiment wherein the wheel speed v n is obtained for each one of the teeth 1 6, the number "j" or the number 
of addresses of the "wheel speed variation" memory is equal to the number of the teeth 16. The address at which the 
wheel speed variation (V - v n ) is stored corresponds to the current number "n". It is also noted that the content at each 
address represents a cumulative wheel speed variation L(V - v n ) = ((V - v n ) at i « 1) + ((V - v n ) at i = 2) + ....) . 

Step S4 is followed by step S5 in which the integer "n" is incremented to repeat steps S3 and S4 if an affirmative 

45 decision (YES) is obtained in the following step S6, that is, if the incremented integer "n" is larger than or equal to "j". If 
a negative decision (NO) is obtained in step S6, this means that the tired wheel 14 has rotated one full revolution, and 
that the wheel speed values v n corresponding to all the teeth 16 have been obtained. With the negative decision (NO) 
obtained in step S6, the control flow goes to step S7 to reset the integer "n" to "1" and increment an integer V. Then, 
step S8 is implemented to determine whether the integer V is smaller than "ISP. If an affirmative decision (YES) is 

so obtained in step S8, the control flow goes back to step S2 to determine again the vehicle speed V, and repeatedly 
implement steps S3 and S4 for obtaining the cumulative wheel speed variations corresponding to all the teeth 1 6 during 
the next one revolution of the tired wheel 14. The wheel speed variation (V - v n ) obtained by each implementation of 
step S4 is added to the content of the address of the "wheel speed variation" memory which corresponds to the integer 
"n". Thus, the content (representative of the cumulative wheel speed variation) of each address of the memory is 

55 increased each time the tired wheel 1 4 is rotated by one revolution. 

In the illustrated embodiment, the vehicle speed V is not updated each time the integer "n" is incremented or each 
time steps S3 and S4 are repeated, on the assumption that the vehicle speed does not change during one full revolution 
of the wheel 14. However, the routine of Fig. 5 may be modified such that the control flow goes back to step S2 each 
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time the affirmative decision (YES) is obtained in step SB, namely, each time the wheel speed variation (V - v n ) corre- 
sponding to each tooth 16 is obtained. 

A negative decision (NO) is obtained in step S8 when the integer Y has become equal to the predetermined number 
"N", namely, when the tired wheel 14 has rotated the predetermined number "N" of revolution. In this case, step S9 is 

5 implemented to divide the content (cumulative wheel speed variation) of each address of the "wheel speed variation" 
memory of the RAM 44, by the number "N", to thereby obtain an average cumulative wheel speed variation Av n . This 
value Av n is an average cumulative variation of the wheel speed v n at each of the teeth 16 of the rotor 10, with respect 
to a reference position of the rotor 10, which is the position of the tooth 16 at which the routine of Fig. 5 is initiated or at 
which steps S3 and S4 are implemented for the first time. Where the routine of Fig. 5 is executed from time to time or 

w continuously during running of the vehicle, the wheel speed calculating/compensating routine of Fig. 6 does not require 
the reference position of the rotor 10 (wheel 14) to compensate the wheel speed for the wheel speed variation. Where 
the routine of Fig. 5 is executed only once immediately after the assembling of the tired wheel 14 or each time a power- 
on switch of the vehicle has been turned on, the reference position of the rotor 1 0 should be known. To this end, a suitable 
marking indicative of the reference position is provided on the rotor 1 0, and a detector for sensing the marking is fixedly 

is disposed, to detect the reference position. 

Where the marking indicative of the reference position of the rotor 10 is provided, the routine of Fig. 5 to obtain the 
value of the average wheel speed variation Av n corresponding to the individual teeth 1 6 may be executed beginning with 
the tooth 16 at the reference position. Alternatively, the routine may be initiated at a desired position of the rotor 10, and 
the obtained values Av n are converted to those with respect to the reference position. 

20 The number "N" of revolutions of the rotor 10 (wheel 14) for which the average cumulative wheel speed variation 
values Av n are obtained should be large enough to minimize an influence of the road surface condition on the obtained 
values Av n . In the present invention, the number "N" is determined depending upon the vehicle speed V, more specifically, 
the number "N" is determined so as to increase with the vehicle speed V, as indicated in the graph of Fig. 7. 

In step S9. the obtained average cumulative wheel speed variation Av n for each tooth 16 is divided by the vehicle 

25 speed V, to obtain a cumulative wheel speed variation rate AVr/V. The obtained rate AVp/V for each tooth 16 is stored at 
an appropriate address of a "wheel speed variation rate" memory of the RAM 44. The values of this rate AVp/V corre- 
sponding to the teeth 16 represent a rotating speed variation inherent to the specific wheel 14. which arises from man- 
ufacturing and assembling errors of the rotor 10 and wheel 14. These values Av n /V are used in the routine of Fig. 6 to 
compensate or adjust the detected speed v of the wheel 14. 

30 The wheel speed calculating/compensating routine of Fig. 6 is executed successively during running of the vehicle. 
In the present embodiment, the routine is executed to calculate the wheel speed from an average time duration between 
predetermined rises of the rectangular pulses generated from the waveform shaper 1 8 during a predetermined sampling 
time. 

The routine of Fig. 6 is initiated with step S1 1 to detect the first and last rises of the rectangular pulses during the 
35 sampling time, count the number of the rises during the sampling time, and determine the two teeth 16 of the rotor 10 
which correspond to the detected first and last rises of the pulses. Described in detail, each time a rectangular pulse 
rises, an interruption routine is executed to detect the moment of the rise on the basis of an output of a timer incorporated 
in the computer 20. A counter is provided to count the number of the rises which occur during the sampling time. Another 
counter is provided to count the number of the rises as counted from the rise corresponding to the reference position of 
40 the rotor 10. This counter is reset at the reference position, and the content of the counter identifies the tooth 16 corre- 
sponding to the rise of the rectangular pulse generated last. 

Then, the control flow goes to step S12 to calculate the average speed v of the wheel 14 during the sampling time, 
on the basis of an average time interval between the rises of the adjacent rectangular pulses during the sampling time. 
Step S12 is followed by step S13 to compensate or adjust the wheel speed v, on the basis of the wheel speed 
45 variation rates Av n1 /V and Av n2 /V corresponding to the two teeth 16 which correspond to the first and last rises of the 
rectangular pulses. These variation rates are read from the "wheel speed variation rate" memory of the RAM 44, on the 
basis of the two teeth 16 determined in step S1 1. The compensation of the wheel speed v is effected according to the 
following equation (1): 

50 v = (1 +(Av n1 - A Vn2 )/2V)v (1) 

It will be understood from the above equation (1) that the wheel speed v is compensated by an amount equal to 
1/2 of a difference between the rates Av n1 /V and Av n2 W multiplied by the pre-compensated value v. 

The compensated wheel speed v obtained in step S13 of the routine of Fig. 6 is processed by the pre-processing 
55 filter 64. and then used by the disturbance observer 52 indicated in Fig. 4. However, the disturbance observer 52 uses 
an angular velocity corresponding to the wheel speed v. To this end, the compensated wheel speed v obtained in step 
S13 is converted into the angular velocity co R of the rim 28, with the radius R of the tire 26 taken into account. The 
calculated angular velocity cdr is stored in an "angular velocity" memory of the RAM 44. 
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The disturbance observer 52 is arranged according to the dynamic model of the tired wheel 14 as illustrated in Fig. 
3. T ^re will be described the arrangement of the disturbance observer 52. 

the tired wheel 14 is simulated as the dynamic model of Fig. 3 wherein the rim 28 and the belt 30 having the 
re^, active moments of inertia J R and J B are connected by the torsion spring 32 having the spring constant K, the following 
state equations (2), (3) and (4) are established so as to define a linear dynamic system: 



J R co R ' = -K6 rb + T -I 



J B co B ' = K6 



RB " T d 



'RB 



COr - COr 



(2) 
(3) 
(4) 



where, 

15 co R : angular velocity of the rim 28, 

cor': angular acceleration of the rim 28, 

co B : angular velocity of the belt 30, 

co B ': angular acceleration of the belt 30, 

6rq: torsion angle between the rim 28 and the belt 30, 

20 T 1 : driving or braking torque detected by the torque detector 68, 

T d : disturbance torque (which is derived from rolling resistance of the wheel 14 which irregularly or regularly 

occur due to steps or undulation on the road surface, respectively). 

Actually, however, a damper is present between the rim 28 and the belt 30. Since the influence of the damper is 
25 relatively small, the presence of the damper is ignored in the above equations according to the present embodiment. 

The above state equations (2) through (4) may be converted into the following equation (5) using vectors and matri- 
ces: 
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40 



The motion of the tired wheel 14 when the spring constant K of the torsion spring 32 is changed from K to K + 
AK due to a change in the air pressure of the tire 26 is expressed by the following equation (6): 
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(6) 



55 



The amount of change AK of the spring constant K is equivalent to the amount of a disturbance acting on the tire 
26 in the normal condition, which disturbance is expressed by the last term of the right member of the above equation 
(6). It will be understood that the last term of the right member of the equation (6) includes the amount of change AK, 
which varies with a change in the air pressure of the tire 26. In other words, the amount of change in the air pressure 
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of the tire 26 can be estimated by estimating the disturbance as expressed by the last term of the right member of the 
equation (6). The disturbance observer 52 should be adapted to estimate an overall disturbance w which includes the 
disturbance torque T d (due to a change of the road surface condition) of the wheel 1 4 as well as the disturbance due to 
the change in the air pressure of the tire 26. The overall disturbance w is represented by the following equation (7): 

r P 1 r-AK/j « i 

[w] = I^-1/Jb J T - + [ A K/J B J 0rb • • • (7) 



Theoretically, however, only one of the three elements of the disturbance w can be estimated by the observer 52. 
Therefore, only the second element w 2 is estimated according to the following equation (8), and the motion state of the 
75 tired wheel 14 is expressed by the following state equation (9): 

w 2 = (-1/J B )T d + (AK/J B )9 RB (8) 




Thus, the disturbance observer 52 operates according to the above equation (9), to estimate the disturbance as 
one of variables which reflect the state of the system. To use the disturbance w 2 of the equation (8) as one of the variables, 
the following equation (10) is used to approximate the disturbance to be estimated: 

35 

w 2 =0 (10) 

The approximation according to the above equation (10) means approximation (zeroth order approximation) of a 
continuously varying disturbance by successive values which change in steps as indicated in the graph of Fig. 8. This 
ao approximation is sufficiently acceptable if the disturbance estimating speed of the disturbance observer 52 is high enough 
to follow the changing rate or velocity of the disturbance to be estimated. The following equation (11) represents an 
expanded system which includes the disturbance w 2 as one of the system variables according to the equation (10): 
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(11) 



In the above equation (11 ), [co B e RB wj' is the state which cannot be detected. In the system using the disturbance 
55 observer 52, not only the disturbance w 2 but also the variables 09 and 6 RB that cannot be measured can be estimated. 
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For simplification, the above equation (1 1) is broken down into the following vectors and matrices: 
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The minimum-order observer to estimate the state [x b ] = [co B 0 RB w 2 ] T is expressed by the following equation 

20 (12): 

[z p l = [A 21 ] [x J + [A 22 ] [z p ] + [B 2 ] [u] + [G] {[x a l - ([A n ] [xj + [A 12 ] [z p ] + [B, ][u])} (12) 
= ([A 21 ] - [G] [A 1 1 ] ) [x J + ([A 22 ] - [G] [A 12 ] ) [z p ] + [G] [x a ' ] + ([B 2 ] - [G] [B, ] ) [u] 

25 

where, 

[z p ] : estimated value of [xj , 
[Zpl: rate of change of the estimated value [z p ], 
30 [G]: gain which determines the estimating speed of the disturbance observer 52. 

The above equation (12) is expressed by the block diaphragm of Fig. 9, wherein [I] represents a unit or identity 
matrix, while "s" represents a Laplace operator. 

If an error between the true value [xj and the estimated value [z p ] is expressed as [e] = [x b ] - [z p ]r and if a differ- 
35 entiation of the error [e] is represented by [el the following equation (13) is obtained: 

[el=([A 22 ].[G][A 12 ])[e] (13) 

The above equation (13) indicates the estimating characteristic of the disturbance observer 52, and the eigen 
40 values of the matrix ([A22] - [G] [A 12 ]) are the pole of the disturbance observer 52. Accordingly, the estimating speed of 
the disturbance observer 52 increases as the eigen values deviate from the origin in the left half of the S-plane. The 
observer gain [G] is suitably determined depending upon a desired value of the estimating speed of the disturbance 
observer 52. 

The disturbance observer 52 arranged as described above receives as an input the angular velocity to R which is 
45 calculated from the compensated wheel speed v of the tired wheel 14 and the tire radius R by the rim speed calcula- 
tor/compensator 45 as described above. On the basis of the angular velocity co R , the disturbance observer 52 estimates 
the disturbance w 2 as expressed by the above equation (8) due to the amount of change AK of the spring constant K of 
the torsion spring 32. The thus obtained disturbance w 2 is expressed as an estimated disturbance value w 2p . In addition 
to the value wgp, the disturbance observer 52 is adapted to obtain an estimated value co Bp of the angular velocity co B of 
50 the belt 30 and an estimated value e RBp of the torsion angle 6 RB between the rim 28 and the belt 30, which co B and e RB 
cannot be detected or measured. 

The correlation calculator 56 operates to calculate correlation functions, on the basis of the estimated disturbance 
value w 2p and torsion angle e RBp . Using the calculated correlation functions, the normalizer 58 performs a normalizing 
operation and calculates the amount of change of the spring constant K of the torsion spring 32. 
55 Referring next to the flow chart of Fig. 1 1 . there will be described the manner of obtaining the amount of change of 
the spring constant K of the torsion spring 32. 

The routine of Fig. 1 1 is initiated with step S21 to reset the integer V to "1", and reset correlation functions C(w 2p , 
0rb p ) and C(e RBp , 0 RBp ) to "0". The correlation function C(w 2p , 6 RBp ) is a function of cross-correlation between the 
estimated value w 2p of the disturbance w 2 expressed by the above equation (8) and the estimated torsion angle value 
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e RBp> while the function C(6 RBpt 6 RBp ) is a function of auto-correlation of the estimated torsion angle value e RBp . In other 
words, the contents of "cross-correlation** and "auto-correlation" memories of the RAM 50 are cleared in step S21 . 

Step S22 is then implemented to read the present estimated disturbance value w 2p( Q and the present estimated 
torsion angle value 6 RBp(i) . Step S22 is followed by step S23 to calculate a product of the estimated disturbance value 
5 w 2p (i) and the estimated torsion angle value 6 RBp (j), and add the product to the last value of the cross-correlation function 
C(w 2p , 8 RBp ) to thereby update the cross-correlation function C(w 2p , 6R B p). In the first cycle of execution of the routine 
of Fig. 11 , the last value of the cross-correlation function C(w2p, 6 RBp ) is zero, and the above-indicated product w 2 p^ x 
e RBp(i) is stored in the "cross-correlation** memory. 

The control f low then goes to step S24 to obtain the square of the estimated torsion angle value 6 RBp ^, and add 
io the obtained square to the last value of the auto-correlation function C(0 RBp , Q^ Bp ) to thereby update the auto-correlation 
function C(6 RBp , e RBp ). The updated value is stored in the "auto-correlation" memory. 

Step S25 is then implemented to determine whether the integer V has become equal to or larger than a reference 
value "M" or not. Initially, a negative decision (NO) is obtained in step S25, and step S26 is implemented to increment 
the integer V, to repeat steps S22, S23 and S24. 
is An affirmative decision (YES) is obtained in step S25 when steps S22-S24 have been repeated the predetermined 
number of times "M". Thus, one cycle of execution of the routine of Fig. 1 1 is terminated. 

While the above reference value "M" may be a predetermined or fixed integer value according to the present inven- 
tion, this value "M" is variable in the present embodiment, as described later, to assure improved accuracy with which 
the tire pressure is estimated. 

20 After the cross-correlation function C(w 2p , 6 RBp ) and the auto-correlation function C(6 RBp , 6 RBp ) have been obtained 
by the correlation calculator 56, the normalizer 58 obtains a normalized value L K according to the following equation 
(14), and the obtained normalized value L K is stored in an L K memory of the RAM 50. 

L K = C(w 2p , 8 RBp )/C(6 RBp , 6 RBp ) (14) 

25 

The value L K obtained according to the above equation (14) may be expressed by the following equation (15), 
on the basis of the above equation (8): 

L K = (-1/J B )C 0 + AK/J B (15) 

30 

Since the value Co is represented by CfT^, 0RB P )/C(e RBp , e RBp ) and is independent of a change of the spring 
constant K, the value Co may be obtained when the air pressure of the tire 26 is normal. It is noted that the value COdp, 
6 RBp ) is a function of cross-correlation between the estimated value T dp of the disturbance torque T d and the estimated 
value 6 RBp of the torsion angle 6 RB . 

35 The determining processor 62 is adapted to compare the normalized value L K = C(w 2p> 8 RBp )/C(e RBp , 6 RBp ) 
stored in the I_k memory of the RAM 50, with a predetermined reference value L K0 stored in the ROM 49. The reference 
value L K0 is a negative value. If the value L K is smaller than the reference value L K0 , this means that the air pressure of 
the tire 26 is abnormally low, and the display device 66 provides an indication informing the vehicle operator of this fact. 
In this respect, it is noted that the ROM 49 stores a data table representative of a relationship between the value L K and 

40 an amount of change AP of the air pressure P of the tire 26, so that the amount of change AP of the air pressure P can 
be determined on the basis of the value L K and according to the L K -AP relationship. 

It will be understood that the correlation calculator 56, normalizer 58 and determining processor 62 constitute means 
for estimating an amount of change from a basic value indicative of a certain kind of wheel information, such as an air 
pressure of the tire 26, to a current value indicative of the current state of the wheel information. 

45 While respective functions of the disturbance observer 52, correlation calculator 56, normalizer 58 and determining 
processor 62 have been separately described, there will be hereinafter described an operation of these elements 52, 
56, 58 and 62 of the computer 47, which cooperate with each other to execute a routine as illustrated in the flow chart 
of Fig. 10. 

The routine of Fig. 10 is initiated with step S31 to initialize the above-indicated reference value "M", that is, to set 
so the value "M" to a predetermined value stored in the ROM 49. Step S32 is then implemented to reset the integer V to 
"1". In the next step S33, the cross-correlation function C(w 2p , 6 RBp ) between the estimated value w 2p of the disturbance 
w 2 expressed by the above equation (8) and the estimated torsion angle value 6 RBp , and the auto-correlation function 
C(8 RB p. e RBp) of tne estimated torsion angle value 0 RBp are both reset to "0". In other words, the contents of "cross- 
correlation" and "auto-correlation" memories of the RAM 50 are cleared in step S33. It will be understood that steps S32 
55 and S33 correspond to step S21 of the flow chart of Fig. 1 1 . 

The control flow then goes to step S34 to read the wheel speed v obtained from the computer 20 through the pre- 
processing filter 64. 

There will be hereinafter described the operation of the pre-processing filter 64. 
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The pre-processing filter 64 receives an input signal indicative of the wheel speed v, from the rim speed calcula- 
tor/compensator 45 of the computer 20. The signal indicative of the wheel speed v, which has not been processed by 
the pie-processing filter 64, generally has a frequency characteristic as indicated in the graph of Fig. 12. To enable the 
distu/oance observer 52 to estimate the disturbance w 2 and the torsion angle e RB with sufficiently high accuracy, using 
5 the simple system of the tired wheel 1 4 as illustrated in Fig. 3, the wheel speed signal to be received by the disturbance 
observer 52 must be substantially held in a predetermined frequency range. To this end, in the present embodiment, 
the pre-processing filter 64 is provided between the rim speed calculator/compensator 45 and the disturbance observer 
52. 

The pre-processing filter 64 has a filtering characteristic as indicated in the graph of Fig. 14, in which the gain (a 

10 unit indicating the strength of the signal) of the signal which is in the predetermined frequency range (between f 1 and 
f2) is equal to 0 t while the gain of the signal outside the predetermined frequency range becomes negative. Through 
this pre-processing filter 64, only a component of the wheel speed signal which is in the predetermined frequency range 
is transmitted from the rim speed calculator/compensator 45 to the disturbance observer 52. If the wheel speed signal 
received by the pre-processing filter 64 has a frequency characteristic over three frequency ranges Rf i, Rf 2 and Rf 3 , as 

15 indicated in the graph of Fig. 12, only a component of the signal which is within the middle frequency range Rf 2 is 
extracted or picked up. and converted into an output signal having a frequency characteristic as indicated in the graph 
of Fig. 13. This output signal generated by the pre-processing filter 64 is suitably used by the disturbance observer 52 
for estimating the disturbance w 2 and the torsion angle e RB . 

Once the wheel speed v is read from the pre-processing filter 64. the control flow goes to step S35 in the flow chart 

so of Fig. 10. in which step the disturbance observer 52 is adapted to estimate the disturbance w 2 and the torsion angle 
6^. Step S36 is then implemented to calculate the cross-correlation function C(W 2p , e RBp ) and the auto-correlation 
function C(e RBp , e RBp ), respectively. Namely, this step S36 corresponds to step S22 through step S24 in the flow chart 
of Fig. 1 1 . Step S37 is then implemented to determine whether the integer V is equal to or larger than the reference 
value "M" or not. In the first cycle of execution of the routine of Fig. 10, the integer V is equal to "1", and a negative 

25 decision (NO) is obtained in step S37. Then, the integer V is incremented in the next step S38. and the control flow 
goes back to step S34. 

If the integer V becomes equal to or larger than the reference value "M" after repeatedly implementing steps S34- 
S36, an affirmative decision (YES) is obtained in step S37. and step S39 is then implemented to calculate the normalized 
value or correlation ratio L K which is obtained by dividing the cross-correlation function C(w2p, 6 RBp ) by the auto-corre- 

30 lation function C(e RBp . e RBp ). In the next step S40, a variation 6 of the correlation ratio L K is calculated. 

In the normal situation, the correlation ratio L K is supposed to vary by small degrees or to be kept at a constant 
value. However, this ratio 1_k tends to vary to a large extent, in other words, the variation 5 tends to be large when the 
initial or preset reference value "M" is not sufficiently large in view of the current running condition of the vehicle. As 
indicated above, the reference value "rvT indicates the number of the wheel speed values v to be used to calculate the 

35 cross-correlation and auto-correlation functions C(w 2p , 8 RBp ) and C(e RBp , 6 RBp ). Since the variation 5 is decreased with 
an increase in the reference value "M", the reference value "M" is suitably determined in this embodiment, depending 
upon the variation 5 obtained in step S40, so that the determined value "rvT can be used for obtaining the correlation 
ratio L K in the next cycle of execution of the routine of Fig. 10. 

Described in more detail, ROM 49 stores a predetermined relationship between the variation 5 and the reference 

40 value "M" as represented by the graph of Fig. 15, in which the reference value "M" increases by steps from its initial 
value or nominal value, with an increase in the variation 6 of the correlation ratio L K . It will be understood from the graph 
of Fig. 15 that the variation 5 is fixed or does not change once the value 6 exceeds a predetermined maximum value 
6 max- The reference value "M" to be used upon the next execution of step S37 is determined according to this relationship 
indicated in the graph of Fig. 15. 

45 The present embodiment is also adapted such that the judgement as to whether the tire pressure is abnormally low 
or not is inhibited when the variation 6 exceeds the above-indicated maximum value 5 MAX - In other words, the judgement 
as described above is allowed to be made only when the variation 5 is equal to or smaller than the maximum value 6 MAX . 
This is because the correlation ratio L K is not reliable when the variation 6 exceeds the maximum value 5max. and the 
judgement on the tire pressure may involve an error if it is made on the basis of the unreliable value L K . The variation 6 

so may be calculated as a difference between the correlation ratio L K obtained in the current cycle of execution of the routine 
of Fig. 1 0 and thai obtained in the previous cycle. It will be understood that the variation 6 is one example of an amount 
on the basis of which the number of values, such as correlation ratios, used for estimating a desired wheel information, 
such as a tire air pressure, is determined. 

After the variation 6 of the correlation ratio L K is calculated in step S40, step S41 is implemented to determine the 

55 reference value "M" to be used in the next cycle, depending upon the variation 5. Step S42 is then implemented to 
determine whether or not the variation 8 exceeds the predetermined maximum value 8 MAX . If a negative decision (NO) 
is obtained in this cycle, the control flow goes to step S43 to determine whether the correlation ratio L K is smaller than 
the negative reference value L K o or not. If a negative decision (NO) is obtained in this cycle, that is, if the correlation 
ratio L K is not smaller than the value L K0 . step S44 is implemented to determine that the air pressure P of the tire 26 is 
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normal. In this case, a command is generated in step S46 for inhibiting the display device 66 from providing an indication 
relating to the air pressure of the tire 26. The control flow then goes back to step S32. If an affirmative decision (YES) 
is obtained in step S43, that is. if the correlation ratio L K is smaller than the reference value L K0 , step S45 is implemented 
to determine that the air pressure P of the tire 26 is abnormally low. In this case, a command is generated in step S46 

5 for activating the display device 66 to provide an indication informing the vehicle operator of an abnormally low level of 
the air pressure P. The control flow then goes back to step S32. 

If the variation 6 exceeds the maximum value &max* an affirmative decision (YES) is obtained in step S42 and the 
control flow goes back to step S32. In this case, steps S43-46 are skipped, and the determination on the abnormality 
of the air pressure P is not made. 

w It will be understood from the description of the present embodiment that the rotor 10, magnetic pickup 12 and 
waveform shaper 18 cooperates with a portion of the rim speed calculator/compensator 45 assigned to calculate the 
wheel speed v to constitute a detecting device, while the disturbance observer 52 and the parameter calculator 53 
constitute estimating means for estimating the air pressure of the tired wheel 14 as wheel information. It will also be 
understood that the pre-processing filter 64 serves as frequency characteristic modifying means for providing an optimum 

is frequency characteristic of the wheel speed signal as a wheel motion state or parameter to be supplied to the estimating 
means. 

Referring next to Figs. 16-19. there will be described a second embodiment of the present invention in the form of 
a tire pressure alarming apparatus. In these figures, the same reference numerals as used in the first embodiment are 
used for identifying structurally and/or functionally corresponding elements, of which detailed explanation will not be 
20 provided. 

In the present embodiment, the computers 20 and 47 are connected to each other, through the pre-processing filter 
80 and a FFT (high-speed Fourier transform) analyzer 82, as shown in Fig. 1 7. Like the first embodiment, the computer 
20 includes the rim speed calculator/compensator 45, as shown in Fig. 16. The computer 47 stores in the ROM 49 
various control programs including a tire pressure alarming routine as illustrated in Fig. 18, and thus constitute a reso- 
25 nance-point detector 84, air pressure calculator 86 and determining processor 88, as shown in Fig. 16. 

The computer 47 is adapted to estimate the air pressure P of the tired wheel 26 from resonance frequencies, on 
the basis of the fact that the resonance frequency of the wheel speed signal which is in a predetermined frequency 
range is reduced with a decrease in the tire pressure P. 

Since the computer 47 estimates the air pressure P on the basis of the frequency characteristic of the wheel speed 
30 signal, as described above, the FFT analyser 82 is provided between the computers 20 and 47, such that the frequency • 
characteristic of the wheel speed signal is analyzed by the FFT analyzer 82 and then supplied to the computer 47. The 
wheel speed signal may be transmitted from the rim speed calculator/compensator 45 directly to the FFT analyzer 82 
for analysis of the frequency characteristic thereof. In this case, however, a component of the signal which is outside 
the predetermined frequency range is also supplied to the FFT analyzer 82, resulting in reduced analyzing accuracy. In 
35 the present embodiment, therefore, the pre-processing filter 80 is provided between the computer 20 and the FFT ana- 
lyzer 82, such that only a component of the wheel speed signal which is in the predetermined frequency range is trans- 
mitted from the computer 20 to the FFT analyzer 82. through the pre-processing filter 80. 

The computer 47 operates to execute the tire pressure alarming routine as illustrated in Fig. 18, in the manner as 
described below. 

40 Initially, step S5 1 is implemented to read the frequency characteristic of the wheel speed signal from the FFT analyzer 
82. Step S51 is followed by step S52 to detect the resonance frequency f 0 of the signal that is in the predetermined 
frequency range, on the basis of the frequency characteristic read in step S51 . Namely, a portion of the computer 47 
assigned to execute step S52 constitutes the resonance-point detector 84 as shown in Fig. 16. The control flow then 
goes to step S53 to determine whether the detected resonance frequency is smaller than a threshold value f th or not. 

45 Although this determination is similar to that made in step S56 (which will be described later) for determining whether 
the air pressure P is smaller than a reference value Po or not, the reason why these two steps are provided will be 
described later. 

The resonance frequency fo detected in step S52 is referred to as a provisional value, to distinguish it from a final 
value which will be obtained in later steps. If the provisional value of the resonance frequency f 0 is equal to or larger 

so than the threshold value f th in this cycle of execution of the routine, a negative decision (NO) is obtained in step S53, 
and step S54 is then implemented to determine that the detected resonance frequency fo is a final value. Step S54 is 
followed by step S55 to determine the air pressure P of the tired wheel 1 4, depending upon the final value of resonance 
frequency f 0 . That is, the air pressure P is determined on the basis of the determined resonance frequency f 0 , according 
to the relationship between the resonance frequency fo and the air pressure P. which is stored in the ROM 49 of the 

55 computer 47. Step S56 is then implemented to determine whether the determined air pressure P is smaller than the 
reference value Po or not. In this cycle of execution of the routine of Fig. 18. the final value of the resonance frequency 
fo is derived from the provisional value, and the provisional value is determined to be equal to or larger than the threshold 
value f lh in step S53. This means that the air pressure P is equal to or larger than the reference value P 0 . Accordingly, 
a negative decision (NO) is obtained this time in step S56, and it is determined in the next step S57 that the air pressure 
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P Is normal. Step S59 is then implemented to inform the vehicle operator of the normality of the air pressure P by means 
of the display device 66. Then, the control flow goes back to step S51 . 

If an affirmative decision (YES) is obtained in step S53, that is, if the provisional value of the resonance frequency 
f o detected in the current cycle is smaller than the threshold value f th , steps S60-S63 are executed for the following reason. 

s When the provisional detected value of the resonance frequency fo is smaller than the threshold value 1 th , there is 

a possibility that the air pressure P is abnormally low. However, the provisional value of the resonance frequency f o may 
suddenly or unexpectedly change even while the air pressure P is normal. In this case, if the tire pressure P is determined 
to be abnormally low and the vehicle operator is informed of the abnormality each time the provisional value fo is deter- 
mined to be smaller than the threshold value f th , an error may arise in informing the abnormality of the air pressure P, 

10 thus reducing the reliability of the computer 47. In the present embodiment, therefore, the air pressure P is not determined 
to be abnormally low immediately after the provisional resonance frequency value fo is determined to be smaller than 
the threshold value f th . Rather, the final value of the resonance frequency fo is obtained by averaging the last several 
values of the resonance frequency fo stored in the RAM 50. which include the detected value obtained in step S52 of 
the current cycle. The number of these values thus averaged is referred to as a reference value "N". The determination 

15 as to whether or not the air pressure P is abnormally low is finally made in step S56, based on the thus obtained final 
value of the resonance frequency. 

While the above reference value "N" may be a fixed or predetermined value, the value "NT may be varied in the 
present embodiment, depending upon a variation 6 of the provisional detected value of the resonance frequency fo- 
More specifically described referring to the graph of Fig. 19, the reference value "N" is increased with an increase in the 

20 variation 6 until it reaches its maximum value 8 MAX , and is fixed after the variation 6 reaches the maximum value 6 MAX . 
In this embodiment, the variation 6 is another example of an amount on the basis of which the number of values, such 
as resonance frequency values, used for estimating a desired wheel information, such as a tire air pressure, is deter- 
mined. 

It will be understood from the above description that step S53 is a step for provisionally determining whether the air 
25 pressure P is abnormally low, on the basis of the provisional detected value of the resonance frequency fo- On the other 
hand, step S56 is a step for finally determining whether the air pressure P is abnormally low, on the basis of the final 
value of the resonance frequency fo- Thus, these steps S53 and S56 are seemingly similar to each other, but are in 
substance different from each other. 

While the variation 6 exceeds the maximum value 5 MAX in the present embodiment, the provisional value of the 
30 resonance frequency f 0 is obtained in step S53, but the determination of the final value of the resonance frequency f 0 
(S63) is inhibited, and the determination as to the abnormality of the air pressure P (S56) is accordingly inhibited. While 
the variation 6 that is larger than the maximum value 6 MAX is interpreted to indicate a sudden change of the provisional 
value of the resonance frequency fo, it is generally known that the resonance frequency f 0 does not suddenly change 
with a change of the air pressure P of the tired wheel 14. Therefore, the final value of the resonance frequency fo is not 
35 obtained and the abnormality of the air pressure is not determined, when the variation 6 exceeds the maximum value 

6 max- , , , 

Step S60 through step S63 of the routine of Fig. 18 are implemented in the manner as described below. 
Initially, step S60 is implemented to calculate the variation 6 based on the provisional detected values of the reso- 
nance frequency f 0 . More specifically, the variation 6 is calculated as an amount of variation of the provisional value 

40 obtained in step S52 in the current cycle of execution of the routine, from the provisional value obtained in the last cycle. 
Step S60 is followed by step S61 to determine whether the calculated variation 6 exceeds the maximum value 6 MAX or 
not. If the variation 5 does not exceed the maximum value 5 MAX , a negative decision (NO) is obtained in step S61 , and 
step S62 is then implemented to determine the reference value "N" on the basis of the calculated variation 5. Namely, 
the reference value "N" corresponding to the current variation 6 is determined according to the relationship between the 

45 variation 5 and the reference value N, as indicated in the graph of Fig. 19. which relationship is stored in the ROM 49 
of the computer 47. The control flow then goes to step S63 to read the last several values (-N) from a plurality of 
provisional values of the resonance frequency fo that have been stored in the RAM 50. and obtain an average value by 
dividing the sum of these values by the reference value "N". This average value is determined to be the current final 
value of the resonance frequency. The control flow then proceeds to step S55 and the following steps. If the currently 

so obtained variation 5 exceeds the maximum value 5 M ax- an affirmative decision (YES) is obtained in step S61 , and the 
control flow goes back to step S51 . Thus, the estimation of the air pressure P is inhibited in this cycle of execution of 
the routine of Fig. 18 

It will be understood from the above description of the present embodiment that the rotor 10, magnetic pickup 12, 
waveform shaper 18, and a portion of the rim speed calculator/compensator 45 assigned to calculate the wheel speed 
55 constitute a detecting device, and that the FFT analyzer 82, resonance-point detector 84 and air pressure calculator 86 
constitute estimating means for estimating the air pressure on the basis of the frequency characteristic of the wheel 
speed signal. It will also be understood that the pre-processing filter 80 constitutes frequency characteristic modifying 
means for providing an optimum frequency characteristic of the wheel speed signal as a wheel motion state supplied 
to the estimating means. 
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Referring next to Figs. 20 and 21, there will be described a third embodiment of the present invention, which is 
similar to the embodiment shown in Fig. 16, except a portion of the computer 47 assigned to execute a tire pressure 
alarming routine which will be described. In the following description, the same reference numerals as used in the second 
embodiment are used for identifying corresponding elements, of which detailed explanation is not provided. 

5 The computer 47 of the present embodiment operates to execute the tire pressure alarming routine as illustrated 

in Fig. 20, which is different in the following aspects from the control routine of Fig. 18 of the second embodiment. 

The present embodiment is different from the second embodiment in that the final value of the resonance frequency 
fo is always obtained by averaging a predetermined number (= M) of the last provisional values of the resonance fre- 
quency fo obtained in the current and last several cycles of execution of this routine. In this case, the number M is a fixed 

10 or predetermined value. Even with the number M being a fixed value, the present embodiment takes the following meas- 
ure, so as to avoid an error in informing the abnormality of the air pressure P, which error may be otherwise incurred by 
a sudden change of the detected value of the resonance frequency fo. 

Namely, the air pressure P is determined to be norma) each time the final value of the reference frequency fo is 
equal to or larger than the threshold value f th . However, the air pressure P is not always determined to be abnormally 

is low when the final resonance frequency value is smaller than the threshold value f m . After the computer 47 determines 
for the first time that the final resonance frequency value is smaller than the threshold value f th , a number (=N) of final 
values of the resonance frequency are obtained, and the air pressure P is determined to be abnormally low only when 
each and every one of the thus obtained final resonance frequency values is determined to be smaller than the threshold 
value f t h The number "1ST of the final values to be obtained is determined depending upon an amount of change Af 0 that 

20 is a difference between the final value obtained in the current cycle of execution of the routine, and the final value obtained 
last time. That is, when the final value becomes smaller than the threshold value f^, it is assumed that the amount of 
change Af 0 is due to a sudden or unexpected change of the resonance frequency f 0 , and the computer 47 continues to 
determine that the air pressure P is normal, for a period of time in which the sudden change is expected to be eliminated. 
If the final values obtained in this period of time are continuously determined to be smaller than the threshold value f th . 

25 the computer 47 determines that the amount of change Af 0 is truly derived from a change of the air pressure P In this 
manner, a false alarming due to the sudden change of the detected resonance frequency fo can be avoided. It is noted 
that the number "N" is increased in steps with an increase in the amount of change Afo. as indicated in the graph of Fig : 21 . 

As in the first embodiment, the present embodiment is adapted such that the wheel speed variation obtaining routine 
of Fig. 5 is executed each time a power-on switch of the vehicle has been turned on. During a certain period of time 

30 needed to finish execution of this routine, however, the resonance frequency fo cannot be detected with high accuracy, 
and the determination as to whether the air pressure is abnormally low is unreliable. Accordingly, the determination on 
the air pressure is inhibited during the execution of the wheel speed variation obtaining routine, even though the reso- 
nance frequency fo is detected. 

Even after the periodic wheel speed variation is obtained, the determination on the air pressure cannot be made 

35 until the predetermined number (=M) of provisional values of the resonance frequency are obtained to determine the 
first final value. However, the vehicle operator is normally ignorant that the tire pressure alarming device is currently in 
this state of operation and thus not able to accurately determine whether the tire pressure is abnormally low. In the 
present embodiment, therefore, the vehicle operator is informed through the display device 66 that the alarming device 
is in the process of obtaining the periodic wheel speed variation, that is. the device is now in its "LEARNING" mode, or 

40 that the alarming device is in the process of obtaining the first final value of the resonance frequency fo, that is, the 
device is now in its "DETECTING" mode. 

The tire pressure alarming routine will be hereinafter described in detail, referring to the flow chart of Fig. 20. 
The routine is started with step S71 to read a "LEARNING" flag from the computer 20. The "LEARNING" flag is in 
its "ON" state while the periodic wheel speed variation is being obtained, that is, while the routine of Fig. 5 is being 

45 executed, and the same flag is in its "OFF state while the wheel speed variation is not obtained. Step S71 is followed 
by step S72 to read the frequency characteristic of the wheel speed signal which is analyzed by the FFT analyzer 82. 
Step S73 is then implemented to detect the current resonance frequency fo of the wheel speed signal which is in the 
predetermined frequency range, on the basis of the frequency characteristic read in step S72. The detected resonance 
frequency f 0 is stored in RAM 50, and is used in the following steps as a provisional value. The control flow then goes 

so to step S74 to determine whether the number of provisional values which have been obtained is equal to "M" or not. If 
the number "M" of the provisional values have not been obtained, a negative decision "NO" is obtained in step S74, and 
step S75 is implemented to place a "DETECTING" flag provided in the RAM 50 in its "ON" state. The "DETECTING" 
flag is placed in the "ON" state while the number of the provisional values which have been obtained has not reached 
the predetermined value "M", and is placed in its "OFF" state when the number of the obtained provisional values has 

55 reached "M". 

The control flow then goes to step S76 to determine whether the "LEARNING" flag is in the "ON" state or not. If the 
flag is in the "ON" state, an affirmative decision (YES) is obtained in step S76, and the control flow goes to step S77 to 
activate the display device 66 to inform the vehicle operator that the periodic wheel speed variation is being obtained, 
that is, the routine of Fig. 5 is being executed. Thus, one cycle of execution of the routine of Fig. 20 is terminated. If the 
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"LEARNING" flag is in the "OFF" state, on the other hand, a negative decision (NO) is obtained in step S76. and step 
S78is then implemented to determine whether the "DETECTING"flagisin the "ON" stateor not. Since the "DETECTING" 
flag is in the "ON" state in this cycle of execution of the routine, an affirmative decision (YES) is obtained in step S78, 
and step S79 is implemented to activate the display device 66 to inform the vehicle operator that the first final value of 
the resonance frequency fo is being detected. 

An affirmative decision (YES) is obtained in step S74 after repeatedly executing steps S71-S73 until the predeter- 
mined number (=M) of the provisional values of the resonance frequency f 0 are detected. Step S80 is then implemented 
to calculate the first final value of the resonance frequency f 0 by averaging the provisional values obtained in the first to 
M-th cycles of execution of the routine of Fig. 20. Step S80 is followedby step S81 to determine whether the thus obtained 
final value is smaller than the threshold value f th or not. If the final value is equal to or larger than the threshold value 
fth a negative decision (NO) is obtained in step S81. and step S8la is implemented to place an "ABNORMAL" flag in 
its' "OFF" state The "ABNORMAL" flag provided in the RAM 50 is in the "OFF" state to indicate that the air pressure P 
is normal and is in the "ON" state to indicate that the air pressure P is abnormally low. Step S81a is followed by step 
S81b to set the integer "n" to 0. The integer "n" represents the number of the final values which have been obtained 
since a certain point of time, in other words, the number of times an affirmative decision (YES) is obtained in step S81 
consecutively before the integer "n" reaches "N", as described later. Step S84 is then implemented to place the "DETECT- 
ING" flag in the "OFF" state since the final value is obtained in step S80 in this cycle of execution of the routine. Step 
S84 is followed by step S76 and the following steps. 

Since the "LEARNING" flag and the "DETECTING" flag are both in the "OFF" state in the current cycle, negative 
decisions (NO) are obtained in steps S76 and S78. Step S85 is then implemented to determine whether the "ABNOR- 
MAL" flag is in the "ON" state or not Since the "ABNORMAL" flag is in the "OFF" state in this cycle, a negative decision 
(NO) is obtained in step S85, and the control flow goes to step S86 to inform the vehicle operator via the display device 
66 that the air pressure P is currently normal. Thus, one cycle of execution of the routine of Fig. 20 is terminated. 

If the final value of the resonance frequency f 0 obtained in step S80 is smaller than the threshold value f an 
affirmative decision (YES) is obtained in step S81 , and step S88 is then implemented to determine whether the integer 
"n" is equal to 0 or not. Since the integer "n" is 0 when step S88 is implemented for the first time, an affirmative decision 
(YES) is obtained in step S88. and step S89 is then implemented to calculate an amount of change Af 0 of the current 
resonance frequency f 0 from the last obtained value, and determine the reference value "N" on the basis of the calculated 
amount of change Af 0 , as described above. Step S90 is then implemented to determine whether the integer "n" is larger 
than the reference value "N" or not. If a negative decision (NO) is obtained in step S90, step S82 is then implemented 
to place the "ABNORMAL" flag in the "OFF" state. Then, the integer "n" is incremented in step S83. The control flow 
then goes to step S84 and the following steps, and a negative decision (NO) is obtained in step S85 since the "ABNOR- 
MAL" flag is currently placed in the "OFF" state. Then, step S86 is implemented to activate the display device 66 to 
inform the vehicle operator that the air pressure P is currently normal. 

If the integer "n" becomes larger than the reference value "N" after repeatedly executing the routine of Fig. 20, an 
affirmative decision (YES) is obtained in step S90. and step S92 is then implemented to place the "ABNORMAL" flag 
in the "ON" state. In the next step S92, the integer "n" is set to 0 for the next cycle of execution of this routine. The control 
flow then goes to step S84 and the following steps. Since an affirmative decision (YES) is obtained in step S85 this time, 
step S93 is implemented to activate the display device 66 to inform the vehicle operator that the air pressure P is now 

abnormally low. . . ^ ^ 

In the present embodiment, the resonance frequency f 0 is detected irrespective of whether it is difficult for the mag- 
netic pickup 12 to detect the wheel speed v with sufficiently high accuracy, in the cases where the tired wheels 14 are 
about to stop rotating, or where the vehicle is running on a rough road, such as a graveled road, for example. However, 
the detection of the resonance frequency f 0 . and the determination on the tire pressure P may be inhibited under such 
conditions. In this case, the tire pressure alarming device will not inform the vehicle operator of the abnormality of the 
air pressure even if the air pressure is actually abnormally low. To avoid misunderstanding of the vehicle operator, it is 
desirable to indicate on the display device 66 that the vehicle is now in the "PREPARING" mode, to inform the vehicle 
operator that the determination on the air pressure P is currently inhibited since the resonance frequency f 0 cannot be 
accurately detected. 

It will be understood from the above description of the present embodiment that the rotor 10, magnetic pickup 1 2, 
waveform shaper 18, and a portion of the rim speed calculator/compensator 45 assigned to calculate the wheel speed 
v constitute a detecting device, and the FFT analyzer 82, resonance-point detector 84 and air pressure calculator 86 
constitute estimating means for estimating the air pressure on the basis of the frequency characteristic of the wheel 
speed signal. It will also be understood that the pre-processing filter 80 constitutes frequency characteristic modifying 
means for providing an optimum frequency characteristic of the wheel speed signal as a wheel motion state supplied 
to the estimating means. 

Referring next to Figs. 22-26, there will be described a fourth embodiment of the present invention in the form of a 
tire pressure alarming apparatus. As shown in Fig. 22, the present apparatus includes the magnetic pickup 12, waveform 
shaper 18, computer 20. vehicle speed sensor 70 (which is not shown in Fig. 22), pre-processing filter 64, computer 47, 
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torque detector 68 (which is not shown in Fig. 22) and display device 66. In this embodiment, the filtering characteristics 
of the pre-processing fitter 64 may be varied. That is. a filter range of the filter 64, i.e., a frequency range of a wheel 
speed signal which passes through the pre-processing filter 64. is variable. Further, the computer 47 includes a filter 
range adjuster 90, in addition to the disturbance observer 52, parameter calculator 53 and the determining processor 62. 

5 The graph of Fig. 23 schematically shows a frequency characteristic of a wheel speed signal generated by the rim 

speed calculator/compensator 45, as observed when the air pressure P of the tired wheel 26 is kept at a normal level, 
and a frequency characteristic of the same signal observed when the air pressure P is lowered from the normal level. 
A frequency component of the wheel speed signal which corresponds to the middle one of three peaks of each frequency 
curve shown in Fig. 23 is suitably used by the disturbance observer 52 for estimating the air pressure. As is understood 

w from the graph of Fig. 23, an optimum frequency range which covers this frequency component varies with the level of 
the air pressure P, and shifts to the lower range as the air pressure P is reduced. Therefore, the optimum frequency 
range, or the filter range of the filter 64, tends to be undesirably broad, if the range is fixed or predetermined, irrespective 
of changes of the air pressure P. As a result, unnecessary frequency components may be picked up by the pre-processing 
filter 64 to pass therethrough, resulting in reduced estimating accuracy of the disturbance observer 52, and reduced 

75 reliability in determining whether the tire pressure is abnormally low. In this embodiment, therefore, the filter range of 
the filter 64 can be precisely changed in accordance with a change in the air pressure P, thus permitting a considerably 
narrow frequency range of the signal to pass through the filter 64. In the present embodiment, the filter range is changed 
depending upon an amount of change AK of the spring constant K, which amount is estimated by the disturbance 
observer 52. 

20 A control program for executing a tire pressure alarming routine as illustrated in Fig. 24 is stored in the ROM 49 of 
the computer 47. This program is adapted to change the filter range or optimum frequency range of the wheel speed 
signal, in the manner as described below. 

The routine of Fig. 24 is started with step S101 to read the wheel speed v from the pre-processing filter 64. Initially, 
the filter range of the pre-processing filter 64 is provisionally set to a range suitable for the case where the air pressure 

25 P of the tired wheel 14 is at a normal level. In the first cycle of execution of the routine of Fig. 24, therefore, the signal 
indicative of the wheel speed v read in step S101 is in the predetermined frequency range corresponding to the provi- 
sional fitter range. 

Step S101 is followed by step S1 02 to activate the disturbance observer 52 to estimate the disturbance w 2 and the 
torsion angle e RB , on the basis of the wheel speed read in step S101. Step S103 is then implemented to estimate an 
30 amount of change AK of the spring constant of the torsion spring 32. This parameter AK may be obtained by calculating 
correlation functions and normalizing, or by a method of least squares, for example. 

The control flow then goes to step S104 to determine whether the current cycle of execution of the routine of Fig. 
24 is the first cycle after the computer 47 is turned on. Since the routine is executed for the first time, an affirmative 
decision (YES) is obtained in step S104, and step S105 is then implemented to change the filter range of the pre- 
ss processing filter 64. More specifically, the filter range is set to an optimum range suitable for the frequency characteristic 
of the wheel speed signal, on the basis of the amount of change AK estimated in step S1 03. For example, the filter range 
is determined according to a data table stored in the ROM 49, which represents a relationship between the amount of 
change AK and the filter range. It will be understood from the above description that a portion of the computer 47 assigned 
to execute steps Si 04 and Si 05 constitutes the filter range adjuster 90, which in turn constitutes means for changing 
40 the optimum frequency range used for providing an optimum frequency characteristic of the wheel speed signal as the 
wheel motion state. The control flow then goes back to step S101 to read again the wheel speed v. 

In the second cycle of execution of the routine of Fig. 24, a negative decision (NO) is obtained in step S104, and 
step S106 is then implemented to determine whether the air pressure P is lower than the reference value Po or not, on 
the basts of the estimated value of the amount of change AP If a negative decision (NO) is obtained in step S1 06, namely, 
45 if the air pressure P is not lower than the value Po. step S107 is implemented to determine that the tire pressure is 
normal. Then, step S108 is implemented to inform the vehicle operator of this fact through the display device 66. If an 
affirmative decision (YES) is obtained in step S106, namely, if the air pressure P is lower than the reference value Po. 
step S1 09 is implemented to determine that the tire pressure is abnormally low. Step S1 08 is then implemented to inform 
the vehicle operator of this fact via the display device 66. In either case, the control flow then goes back to step S101 . 
so It will be understood from the above description that steps S107-S109 for determining the abnormality of the air 

pressure P is inhibited when the instant routine is executed for the first time, and is first allowed in the second cycle. 
This is because the wheel speed v supplied to the disturbance observer 52 is obtained through the pre-processing filter 
64 having the provisional filter range, and the amount of change AK estimated on the basis of this wheel speed v may 
not be sufficiently reliable, and thus cannot be relied upon for determining the abnormality of the air pressure. 
55 In the above-described routine of Fig. 24, the filter range or optimum frequency range is always adjusted in the first 
cycle of execution, while inhibiting the determination as to whether the air pressure is abnormally low. However, the 
determination need not be inhibited even in the first cycle where the air pressure can be estimated with high accuracy. 
Accordingly, the illustrated routine may be modified such that the determination on the air pressure is made without 
changing the filter range even in the first cycle of execution of the routine, where it is determined that the air pressure 
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P can be estimated with high accuracy. The thus modified routine is illustrated in the flow chart of Fig. 25, and will be 
hereinafter described in detail. 

In the present routine, steps S1 1 1 -S1 14 are equivalent to steps S101 -S104 of the previous routine of Fig. 24. If this 
routine is executed for the first time, an affirmative decision (YES) is obtained in step S114. and step S115 is then 
implemented to determine whether an absolute value of the amount of change AK (an amount of deviation from the 
nominal value K) is equal to or smaller than the reference value AKo. When the absolute value of the amount of change 
AK is equal to or smaller than the reference value AKo, the disturbance observer 52, which receives the wheel speed 
signal in the provisional frequency range from the pre-processing filter 64, generally assures sufficiently high estimating 
accuracy If an affirmative decision (YES) is obtained in step S1 15. step S1 17 and the following steps are then imple- 
mented to determine whether the air pressure P is abnormally low. Thus, the determination on the air pressure may be 
made even in the first cycle of execution of the routine. Since the absolute value of the amount of change AK is equal 
to or smaller than the reference value AK 0 as determined in step S1 1 5. an amount of change of the tire pressure P from 
the nominal pressure is accordingly small, and a negative decision (NO) is obtained in step S1 17. Step S118 is then 
implemented to determine that the tire pressure is normal, and this determination is displayed in the next step S1 19. if 
the absolute value of the amount of change AK obtained in this cycle is larger than the reference value AKo. a negative 
decision (NO) is obtained in step S1 1 5. and the control flow then goes to step S1 16 to change the filter range as in step 
S105 of the routine of Fig. 24. while inhibiting in the first cycle the determination as to whether the air pressure P is 

^twHIbe understood that a portion of the computer 47 assigned to execute steps S1 14-S1 16 of this routine constitute 
the filter range adjuster 90. which in turn constitutes means for changing the optimum frequency range used for providing 
an optimum frequency characteristic of the wheel speed signal as the wheel motion state, which is to be supplied to the 

estimating means. , .. . 

The filter range or optimum frequency range is not changed in the second and subsequent cycles of execution of 
each of the above-described two routines of Fig. 24 and Fig. 25. However, the control routine may be modified as illus- 
trated in the flow chart of Fig. 26. such that the filter range is changed in each of the second and subsequent cycles as 

In^pTesenUoutine. steps S121-S124 are equivalent to steps S101-S104 or steps S11 1-S1 14 of the routines of 
Figs 24 and 25 If this routine is executed for the first time, an affirmative decision (YES) is obtained in step S124. and 
step S1 25 is then implemented to change the filter range of the pre-processing filter 64. while inhibiting the determination 
as to whether the air pressure P is abnormally low. _ _ nVKtainBrt m 

In the second and subsequent cycles of execution of the present routine, a negative decision (NO) is obtained in 
step S124 Step S126 is then implemented to calculate a deviation AAK of the estimated amount of change AK obtained 
in the current cycle, from the estimated amount of change AK obtained in the previous cycle, and determine whether an 
absolute value of the deviation AAK is equal to or smaller than a positive reference value AAKo or not. Namely, the 
stability of the estimated value AK obtained in the current cycle is evaluated in this step S126. If an affirmative decision 
(YES) is obtained in step S1 26. that is, if the absolute value of the deviation AAK is equal to or smaller than the reference 
value AAKo this means that the estimated value AK currently obtained is sufficiently stable and reliable, and the current 
filter range 'corresponds to an optimum frequency range suitable for the frequency characteristic of the wheel speed 
signal Then steps S128-S130 are implemented to determine whether the air pressure P is abnormally low. 

If the absolute value of the deviation AAK obtained in this cycle is larger than the reference value AAKo. a negative 
decision (NO) is obtained in step S126, and the filter range is changed in step S127. More specifically, if the deviation 
AAK is a negative value, which indicates that the air pressure P is lowered, the whole filter range is shifted by a prede- 
termined amount At to a lower frequency range. If the deviation AAK is a positive value, which indicates that the air 
pressure P is increased, the whole filter range is shifted by a predetermined amount Af to a higher frequency range. As 
a result of repeatedly executing steps S126. S127 and S121-124, therefore, the absolute value of the deviation AAK 
becomes equal to or smaller than the reference value AAKo, and an affirmative decision (YES) is eventually obtained in 
step S126. At this point of time, the filter range corresponds to an optimum frequency range suitable for the frequency 
characteristic of the wheel speed signal. 

It will be understood from the above description that a portion of the computer 47 assigned to execute steps S126 
and S1 27 constitutes the filter range adjuster 90, which in turn constitutes means for changing the optimum frequency 
range so that an amount of variation of the spring constant K is not greater than a predetermined value. 

Referring next to Figs. 27-35. there will be described a fifth embodiment of the present invention. 

A frequency characteristic of the wheel speed signal generated by the rim speed calculator/compensator 45 is 
represented by a curve having a plurality of peaks, as shown in the graph of Fig. 1 2. The frequency range (Rf 0 to which 
the peak having the lowest frequency belongs involves information relating to an amount of change AJ of the moment 
of inertia J and the frequency range (Rfz) to which the peak having the second lowest frequency belongs involves 
information relating to an amount of change AK of the spring constant K. Accordingly, the accuracy of the disturbance 
observer 52 of estimating the amount of change AJ of the moment of inertia J is improved with an increase in the strength 
of a frequency component that is represented by the peak having the lowest frequency. Similarly, the accuracy of the 
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disturbance observer 52 of estimating the amount of change AK of the spring constant K is improved with an increase 
in the strength of a frequency component that is represented by the peak having the second lowest frequency. 

The tired wheels 14 are connected to a vehicle body through a suspension system, such that the wheels 14 are 
displaceable relative to the vehicle body. Where the motion of each wheel 14 relative to the vehicle body through the 

5 suspension is observed in a transverse direction of the vehicle, the motion locus of the center of the wheel 1 4 is generally 
an arc, rather than a straight line which extends in the vertical direction of the vehicle. 

Fig. 27 shows an example of a double wishbone type rear suspension system, in which a vehicle frame 100 is 
connected to a carrier 102 which rotatably supports each of the wheels 14, through an A-shaped upper arm 104 and a 
pair of lower arms 1 06, such that the wheel 1 4 is displaceable relative to the vehicle frame 1 00. In this case, the center 

10 of the wheel 14 has a motion locus as illustrated in Fig. 28(b), for example. The suspension system shown in Fig. 27 
further includes a shock absorber 1 08, a coil spring 1 1 0, an upper support 112, and a strut bar 116, with respect to each 
wheel 14, and a stabilizer 1 14. 

The motion locus of the point which connects the carrier 102 and the upper arm 104 is supposed to lie on the same 
circle on which the motion locus of the point which connects the carrier 102 and lower arms 106 lies. In this case, the 

75 instantaneous center O, which is the center of the above circle, is positioned apart from the wheel 14 toward the front 
side of the vehicle. As a result, the wheel 14 vibrates relative to the vehicle body not only in the vertical direction but 
also in a front-to-rear direction, in other words, in a running or longitudinal direction of the vehicle. Therefore, the wheel 
1 4 moves relative to a road surface in the front-to-rear direction at a point of contact between the wheel 1 4 and the road 
surface. Since the wheel speed signal, which is hardly influenced by the vertical vibrations of the wheel 14, is strongly 

20 influenced by the vibrations applied in the front-to-rear direction, the front-to-rear vibrations appear in the wheel speed 
signal as vibrations occurring at an unsprung resonance point On the other hand, a frequency range which substantially 
coincides with the unsprung resonance point includes a frequency component of the wheel speed signal as represented 
by the peak having the lowest frequency in the graph of Fig. 12. The strength of this frequency component influences 
the accuracy of estimating the moment of inertia J. In sum, there is a predetermined relationship between the motion 

25 locus of the center of the wheel and the frequency characteristic of the wheel speed signal, in particular, the frequency 
characteristic at the unsprung resonance point. Thus, an optimum or appropriate motion locus can be obtained by utilizing 
this relationship, so as to obtain an optimum frequency characteristic of the wheel speed signal, thus eliminating a need 
to provide the pre-processing filter as used in the preceding embodiments. 

In view of the above, the inventors of the present invention designed three kinds of suspensions having different 

30 geometries, which include: a suspension as shown in Fig. 28B, in which the instantaneous center O is spaced from the 
center of the wheel by a generally known distance; a suspension as shown in Fig. 28A, in which the instantaneous center 
O is shifted toward the front side of the vehicle and located almost indefinitely far from the position shown in Fig. 28B; 
and a suspension as shown in Fig. 28C, in which the instantaneous center O is located nearer to the wheel 14, as 
compared with the position shown in Fig. 28B. 

35 The geometry of the suspension as shown in Fig. 28A may be achieved by positioning the upper arm 104 and the 
lower arms 1 06 relative to each other, as shown in Fig. 29, for example. In this case, the center of the wheel has a motion 
locus as indicated in Fig. 31 by the solid line ® . The geometry of the suspension as shown in Fig. 28B may be achieved 
by positioning the upper arm 104 and the lower arms 106 relative to each other, as shown in Fig. 30, for example. In 
this case, the center of the wheel has a motion locus as indicated in Fig. 31 by the one-dot chain line © . With the 

40 geometry of the suspension as shown in Rg. 28C, the center of the wheel has a motion locus as indicated in Fig. 31 by 
the broken fine ® . 

With each of the above three suspensions installed on a motor vehicle, the present inventors obtained frequency 
characteristics of wheel speed signals generated by the respective vehicles. As is apparent from the results as shown 
in the graph of Fig. 32, the power spectrum (strength) of the frequency component of the wheel speed signal, which 

45 substantially corresponds to the unsprung resonance point, increases as the instantaneous center O is located closer 
to the center of the wheel and the motion locus has a reduced radius of curvature, with a result of increased wheel 
vibrations in the front-to-rear direction. Since the power spectrum at the unsprung resonance point influences the accu- 
racy of estimating the moment of inertia J, as described above, the accuracy for estimating the moment of inertia J is 
improved by locating the instantaneous center O as close as possible to the center of the wheel. 

so On the other hand, the power spectrum (strength) of the frequency component of the wheel speed signal, which 
substantially corresponds to the unsprung resonance point, decreases with increases in the distance between the instan- 
taneous center O and the wheel center, and the radius of curvature of the motion locus, which lead to reduced wheel 
vibrations in the front-to-rear direction. This results in an increased power spectrum of a frequency component located 
next to the above frequency component whose power spectrum is reduced. As described above, the accuracy of esti- 

55 mating the spring constant K is influenced by the power spectrum of the frequency component as represented by the 
peak having the second lowest frequency in the graph of Fig. 12. Thus, the accuracy of estimating the spring constant 
K and eventually the air pressure is improved by locating the instantaneous center O as far as possible from the wheel 
center. 
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Then, the wheel information estimating apparatus of the present embodiment as illustrated in Fig. 33 was installed 
in both the vehicle having the suspension of Fig. 29 with the quasi-straight locus indicated by the solid line © in Fig. 
31 , and the vehicle having the suspension of Fig. 30 with the curved locus indicated by the one-dot chain line © in Fig. 
31. In the vehicle having the suspension with the quasi-straight locus, the disturbance observer 52 operates to estimate 
an amount of change AK of the spring constant, as a wheel information or parameter of the tired wheel 1 4. In the vehicle 
having the suspension with the curved locus, on the other hand, the disturbance observer 52 operates to estimate an 
amount of change AJ of the moment of inertia, as a wheel information or parameter of the wheel 14. 

One example of the disturbance observer 52 adapted to estimate the amount of change AK of the spring constant 
has been described with respect to the first embodiment. There will be hereinafter described another example of the 
disturbance observer 52 adapted to estimate the amount of change AJ of the moment of inertia of the tired wheel 14. 
The moment of inertia of the wheel 14 may be either a moment of inertia J R of the rim 28 or a moment of inertia J B of 
the belt 30. 

There will be first described in detail the manner of estimating an amount of change AJ B of the moment of inertia 
J B of the belt 30, by the disturbance observer 52. 

The motion of the tired wheel 14 when the moment of inertia J B has changed by AJ B is expressed by the following 

equation (16): 
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The disturbance to be estimated is the second element of the last term of the right member of the above equation 
(16). If the disturbance w 2 is defined by the following equation (17), an expanded system as represented by the above 
equation (1 1) is obtained from the following state equation (18): 



w 2 = (-1/J B )T d -(AJ B /J B )co B ' 
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Thus, the disturbance observer 52 operates to estimate the disturbance due to the change AJ B of the moment 
of inertia J B of the belt 30. 

The change AJ B of the moment of inertia of the belt 30 is then obtained by executing a routine for calculating cor- 
relation functions for variation of the moment of inertia J B , as illustrated in the flow chart of Fig. 34. 

In sum, a cross-correlation C(w 2p , cd^*) and an auto-correlation function 0(0^', co Bp ) are obtained from the esti- 
mated value w 2p of the disturbance w 2 expressed by the above equation (1 7) and the estimated angular acceleration 
o) B p of the belt 30. Then, a normalized value L JB is calculated by the following equation (1 9) by the normalizer 58, and 
the calculated normalized value L JB is stored in an L JB memory of the RAM 50. 



-JB =C(W 2p , W Bp yC(cO Bp \ ©Bp") 



(19) 
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The value L JB obtained according to the above equation (19) may be expressed by the following equation (20), 
on the basis of the above equation (17): 

Lj B = (-1/J B )Ci-AJ B /J B (20) 

5 

The value Ci means 0(7^, a^pVCfcDBp', o^'). 
The amount of change AJ B corresponding to the thus calculated value L JB is obtained according to a relationship 
between the value L JB and the amount of change AJ B , which is stored in the ROM 49. 

Then, the manner of estimating the amount of change AJ R of the moment of inertia J R of the rim 28 by the disturbance 
10 observer 52 will be explained. 

Normally, the moment of inertia J R of the rim 28 will not change. If the currently used metal wheel member 24 is 
replaced by a new one. the moment of inertia J R changes from one value to another. If the same moment of inertia J R 
were used even after the new wheel member 24 is used, there would be an error in the estimation of the overall distur- 
bance acting on the tired wheel 14. In the present embodiment, therefore, the disturbance observer 52 is arranged to 
75 also estimate the disturbance due to a change of the moment of inertia J R of the rim 28. 

The motion ef.the tired wheel 14 when the moment of inertia J R has changed by AJ R is expressed by the following 
equation (21): 
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The disturbance to be estimated is the first element of the last term of the right member of the above equation 
(21). If the disturbance is defined by the following equation (22), an expanded system as represented by the following 
equation (24) is obtained from the following state equation (23): 

35 w 1 =-(AJ R /J R )co R ' (22) 
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Thus, the disturbance observer 52 operates to estimate the disturbance due to the change AJ R of the moment 
of inertia J R of the rim 28. 

The change AJ R of the moment of inertia of the rim 28 is then obtained by executing a routine for calculating cor- 
relation functions for variation of the moment of inertia J Rl as illustrated in the flow chart of Fig. 35. 

In sum, a cross-correlation C(w 1p , © R ') and an auto-correlation function C(a> R \ co R ) are obtained from the estimated 
value w 1p of the disturbance Wi expressed by the above equation (22) and the estimated angular acceleration a> R of 
the rim 28. Then, a normalized value L JR is calculated by the following equation (25) by the normalizer 58, and the 
calculated normalized value L JR is stored in an L JR memory of the RAM 50. 



L JR = C(w 1p , co R ')/C(© R \ co R ) 



(25) 



The value L JR obtained according to the above equation (25) may be expressed by the following equation (26), 
40 on the basis of the above equation (22): 



L JR = -AJ R /J R 



(26) 



The amount of change AJ R corresponding to the thus calculated value L JR is obtained according to a relationship 
between the value L JR and the amount of change AJ R , which is stored in the ROM 49. 

When the disturbance observer 52 estimates the moment of inertia J B of the belt 30, the determining processor 62 
of the computer 47 determines whether the estimated amount of change AJ B is larger than a predetermined positive 
reference value AJ B0 . If an affirmative decision (YES) is obtained, a suitable indication is provided on the display device 
66 to inform the vehicle operator that a foreign matter is inserted or stuck in the tire 26, or a chain is attached to the tire 
26. The determining processor 62 also determines whether the estimated amount of change AJ B is smaller than a 
predetermined negative reference value AJ B1 . If an affirmative decision (YES) is obtained, a suitable indication is pro- 
vided on the display device 66 to inform the vehicle operator that the amount of wear of the tire 26 exceeds a tolerable 
upper limit. 

When the disturbance observer 52 estimates the moment of inertia J R of the rim 28. on the other hand, the deter- 
mining processor 62 of the computer 47 determines whether the absolute value of the estimated amount of change AJ R 
is larger than a predetermined positive reference value AJ R0 . If an affirmative decision (YES) is obtained, a suitable 
indication is provided on the display-device 66 to inform the vehicle operator that the wheel metal member 24 was 
replaced by a wrong or inappropriate one. 



26 



BNSDOCID <EP 0699546A1_I_> 



EP0 699 546 A1 

Thus, when the disturbance observer 52 is adapted to estimate the spring constant K and the air pressure P of the 
tired wheel 14, the suspension is designed by suspension adjusting means, such that the upper arm 104 and the lower 
arms 1 06 are disposed parallel to each other as seen in the transverse direction ol the vehicle, whereby the instantaneous 
center O is located further away from the wheel 14 than the nominal position of the center O. 

5 When the disturbance observer 52 is adapted to estimate the moment of inertia J of the tired wheel 1 4, on the other 

hand, the suspension is designed by suspension adjusting means, such that the upper arm 104 and the lower arms 1 06 
form an acute angle as seen in the transverse direction of the vehicle, whereby the instantaneous center O is located 
nearer to the wheel 14 than the nominal position of the center O. 

Referring next to Fig. 36, there will be described a sixth embodiment of the present invention, which includes the 

w magnetic pickup 12, waveform shaper 18, computer 20, vehicle speed sensor 70 (which is not shown in Fig. 36), pre- 
processing filter 64, computer 47 and display device 66, as in the embodiment shown in Fig. 2. However, the respective 
functions of the disturbance observer 52, parameter calculator 53 and determining processor 62 of the computer 47 
shown in Fig. 35 are different from those of the embodiment of Fig. 2, as described below. 

The present embodiment is adapted for a tire model which is a simple body of rotation having the moment of inertia 

15 J, as shown in Fig. 37. This tire model represents a tire-suspension model as well as a suspension model. In Fig. 37, it 
is to be noted that 



T d : disturbance torque received from the road surface 

Tt : driving or braking torque acting on the wheel 1 4 

20 co: angular velocity of the wheel 14 (which is equal to the angular velocity of the rim 28) 

K. : equivalent spring constant of the suspension 

Ds: equivalent damping coefficient of the suspension 



It is generally known that a locus indicative of the motion of the tired wheel 14 connected to the vehicle body 
25 through the suspension, as seen in the transverse direction of the vehicle, includes a component which extends in the 
front-to-rear direction or longitudinal direction of the vehicle, as well as a vertical component. Accordingly, the tire-sus- 
pension model shown in Fig. 37 is designed as a combination of vibrations of the wheel 14 due to the suspension, and 
the rotary motion of the wheel 14 itself. 

The equation of motion of the above model is as follows: 

30 

Jco' = -K s Jcodt-D s co + T 1 -T d (27) 

where, 

35 o>': angular acceleration of the wheel 1 4. 

If only vibrations among various motions of the wheel 1 4 and the suspension are taken into consideration, the driving 
or braking torque T 1 may be considered as a fixed value. Thus, the above equation (27) may be converted into the 
following equation (28): 

40 

Jco' = -K s Jcodt-D s co-T d (28) 

The following equation (29) is then obtained by differentiating both the right and left members of the equation 
(28) by time T: 

45 

Jco" = -K s co-D s co'-T d ' (29) 

where, 

so co" : differentiated value of angular acceleration co* 

Td': differentiated value of disturbance torque. 
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The above equation (29) is converted into the following state equation (30): 
G)' If 0 1 

co" J L-Ks /J -D s /J 



]C-] + [-T°'/J--- (30) 



10 



The spring constant Ks depends on the rigidity of the suspension, and is thus considered not to substantially 
change. It is therefore assumed that the moment of inertia J changes due to a change of the air pressure P. and the 
damping coefficient D s changes depending upon the rate of change of the angular velocity o>, which changes due to a 
change of the moment of inertia J. 

The above equation (30) is converted into the following equation (31 ) when the moment of inertia J and the damping 
15 coefficient D s have changed by AJ and ADs. respectively: 



20 



Leo" J 



0 



-Ks /J -Ds /J 



■ ] t w:i 



• (31) 



where, 



25 



30 



40 



45 



50 



AJKs AJDs - J AD$ , 1_ 

w — - qj + — — - CO 



J (J + AJ) J (J + A J) J +A J 



That is, changes of the moment of inertia J by AJ and the damping coefficient D s by AD S are equivalent to the 
disturbance w applied to the tired wheel 14 in its normal state, which disturbance w is represented by the last term of 
35 the equation (31). The disturbance observer 52 then estimates the disturbance w and the angular acceleration co\ on 
the basis of the angular velocity co as an output signal of the pre-processing filter 64, which consists only of a varying 
component of the angular velocity co as an input signal of the pre-processing filter 64. A delay of the disturbance observer 
52 in estimating the above parameters can be reduced to infinitesimal. If the delay is ignored at present, an estimated 
value of the disturbance w is represented by the following equation (32): 



AJKs AJDs -JADs , 
w — . — — — co 

J (J + AJ) J (J + AJ) 

1 



— T d ' • ■ • (32) 

J+AJ 



The parameter calculator 53 calculates the amount of change AJ of the moment of inertia J, on the basis of the 
above equation (32), using the estimated value of the disturbance w, estimated value of the angular acceleration co', and 
55 detected value of the angular velocity co. in the manner as described below. 
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The above equation (32) can be converted into the following equation (33): 

AJKs -JADs 



W P = 



AJKs 
J (J+AJ) 



J (J + AJ) 



it:.] 



10 
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J+AJ 
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= [e ' e ° ] C}--^7ZT T '' 
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= e T r- 
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J+AJ 



T, ' 



(33) 
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where, 
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J (J + AJ) 
AJKs -JADs 
L J (J + AJ) 



40 



On the basis of the above equation (33), a column vector 0 as an unknown parameter is estimated by a method 
of least squares. The column vector 0 is represented by the following equation (34), using the above definition of 0: 
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Thus, the amount of change AJ of the moment of inertia J is represented by the following equation (35): 

j 2 e,, . . . (35) 

A J = 

Ks - J ©ip 

The rate of change of the moment of inertia J is represented by the following equation (36): 

= • • ■ (36) 

J Ks - je« P 



The determining processor 62 determines whether anything is wrong with the tired wheel 1 4, on the basis of the 
amount of change AJ of the moment of inertia received from the parameter calculator 53. When the determining proc- 
essor 62 determines that the wheel 1 4 is not in the normal state, the display device 66 is activated to inform the vehicle 
operator of this fact. 

The wheel speed signal generated by the rim speed calculator/compensator 45 generally has a frequency charac- 
teristic as shown in the graph of Fig. 38. The disturbance observer 52 estimates the amount of change AJ of the moment 
of inertia J with high accuracy, with respect to an optimum frequency range of the signal which is defined by f A and f 2 in 
the graph of Fig. 38, for example. In the present embodiment, the pre-processing filter 64 is adapted to pick up only a 
component of the wheel speed signal which is in the predetermined optimum frequency range, so that the wheel speed 
signal received by the disturbance observer 52 has a frequency characteristic as represented by the graph of Fig. 39. 
for example. Accordingly, the disturbance observer 52 of the present embodiment can adopt such a simple model for 
the tired wheel 14 as illustrated in Fig. 37, assuring sufficiently high estimating accuracy. 

Referring next to Figs. 41 -43, there will be described a seventh embodiment of the present invention. 

It is recognized that vibrations of the tired wheel 14 have a strong relation with a degree of contact between the 
wheel 14 and the road surface, in other words, a road holding characteristic of the wheel 14. That is, the wheel 14 is 
less likely to vibrate when it is in good contact with the road surface. It is also recognized that vibrations of the wheel 1 4 
generally include those applied in the front-to-rear direction or running direction of the vehicle, as well as those applied 
in the vertical direction, as described above. Therefore, the wheel speed signal indicative of the rotating speed of the 
wheel 1 4 is influenced by the vibrations of the wheel 1 4. More specifically, the frequency component of the wheel speed 
signal which is in a range of the unsprung resonance frequency has a power spectrum (strength) which increases with 
an increase in the unsprung vibrations. The graph of Fig. 41 represents a frequency characteristic of the wheel speed 
signal obtained while the vehicle is running with the wheels 14 being in a bad contact with the road surface. The graph 
of Fig. 42 represents a frequency characteristic of the wheel speed signal obtained while the vehicle is running with the 
wheels 14 being in a good contact with the road surface. 

In view of the above recognitions, the wheel information estimating apparatus of the present embodiment is adapted 
to estimate the degree of contact of the wheel 14 with the road surface, on the basis of the wheel speed signal v. As 
shown in Fig. 40, the estimating apparatus includes the magnetic pickup 12, waveform shaper 18, computer 20, vehicle 
speed sensor 70 (not shown), and pre-processing filter 92, as in the embodiment of Fig. 17. The present apparatus 
further includes a FFT analyzer 94, and a computer 95 which replaces the computer 47 of the embodiment of Fig. 17. 
The computer 95 stores in the ROM 49 various control programs including a routine for determining a degree of contact 
between the wheel 1 4 and the road surface, or a road holding characteristic of the wheel 1 4, as illustrated in the flow 
chart of Fig. 43, and thus constitutes an averaging processor 96 and a determining processor 98 as shown in Fig. 40. 
A suspension control device 120 is connected to the computer 95, and receives therefrom the estimated degree of 
contact between the wheel 14 and the road surface. The suspension control device 120 utilizes the estimated degree 
of contact, to control a particular characteristic of a suspension 124 of the vehicle, by means of an actuator 122. Where 
the device 120 is adapted to control the damping characteristic of the shock absorber as the particular characteristic of 
the suspension 1 24, for example, the actuator 1 22 may take the form of a motor for changing the flow resistance in the 
shock absorber. There will be hereinafter described each of significant elements of the wheel information estimating 
apparatus of Fig. 40. 

The pre-processing filter 92 functions to pick up only a frequency component of the wheel speed signal received 
from the rim speed calculator/compensator 45, which component is in a range of the unsprung resonance frequency. 
The thus selected frequency component is transmitted to the FFT analyzer 94 which calculates the frequency charac- 
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teristic of the wheel speed signal received from the pre-processing filter 92. The thus calculated frequency characteristic 
is transmitted to the computer 95. 

The averaging processor 96 averages a plurality of results of frequency analysis, which have been sequentially 
supplied from the FFT analyzer 94, so as to improve the accuracy of analyzing the frequency characteristic. The deter- 

s mining processor 98 utilizes the result of the frequency analysis finally obtained by the averaging processor 96, and thus 
determines, on the basis of the result, the maximum value G v of gain (which is a unit of the strength of the signal) in the 
frequency range of the unsprung resonance frequency. When the maximum gain Gv is equal to or larger than a prede- 
termined reference value G^q, the determining processor 98 determines that the tired wheel 14 is in bad contact with 
the road surface. When the maximum gain Gy is smaller than the reference value G^, the processor 98 determines that 

w the wheel 14 is in good contact with the road surface. The suspension control device 120 keeps watching the determi- 
nation of the determining processor 98, and makes the damping characteristic of the suspension 124 harder, through 
the actuator 122, when the wheel 1 4 is determined to be in a bad contact with the road surface. This assures improved 
stability of the vehicle, which might be otherwise deteriorated due to such a bad contact of the wheel 14 with the road 
surface. 

is The averaging processor 96 and determining processor 98 operate to execute the routine of determining a degree 
;of contact between the wheel 14 and the road surface, as illustrated in the flow chart of Fig. 43 and described below. 

The routine is started with step S201 to clear the content of a frequency characteristic memory of the RAM 50. The 
frequency characteristic memory stores the result of analysis of the frequency characteristic of the wheel speed signal 
supplied from the FFT analyzer 94. Step S202 is then implemented to set the integer "n n to 1 . The integer "n" represents 

20 the number of frequency characteristics which are currently stored in the frequency characteristic memory. The control 
flow then goes to step S203 to read a frequency characteristic of the wheel speed v supplied in this cycle from the FFT 
analyzer 94, and stores it in the frequency characteristic memory. The integer "n" is incremented in the next step S204, 
and step S205 is then implemented to determine whether the integer "n" is equal to or larger than a predetermined 
reference value "N", i.e. . whether all the frequency characteristics needed to effect the subsequent averaging processing 

25 in the step S206 have been obtained. If the integer "n" is smaller than the reference value "N", a negative decision (NO) 
is obtained in step S205. and the control flow goes back to step S202. 

When the integer "n" becomes equal to or larger than the reference value "N" after repeatedly executing steps S202- 
S204, an affirmative decision (YES) is obtained in step S205, and the averaging processing is effected in step S206. 
Namely, the maximum gain Gv of the signal in the unsprung frequency range is calculated on the basis of the predeter- 

30 mined number of frequency characteristics which have been stored in the frequency characteristic memory. Then, step 
S207 is implemented to determine whether the calculated gain Gv is equal to or larger than the reference value G^. If 
an affirmative decision (YES) is obtained in step S207, step S208 is implemented to place a "BAD CONTACT" flag 
indicative of a bad contact of the wheel 14 with the road surface in its "ON" state, thereby to indicate that the wheel 14 
is in a bad contact with the road surface, if the gain is smaller than the reference value Gv 0 , a negative decision (NO) 

35 is obtained in step S207, and step S209 is implemented to set the "BAD CONTACT flag in its "OFPstate, thereby to 
indicate that the wheel 1 4 is in a good contact with the road surface. In either case, the control flow goes back to step S201 . 

While the FFT analyzer 94 is provided separately from the computer 95 in the present embodiment, the analyzer 
94 may be incorporated in the computer 95. 

In the present embodiment, the FFT analyzer 94 does not receive the wheel speed v directly from the rim speed 

40 calculator/compensator 45. Rather, the pre-processing filter 92 provided between the rim speed calculator/compensator 
45 and the FFT analyzer 94 is adapted to pick up or select only the necessary frequency component of the wheel speed 
v, namely, the component in the unsprung resonance frequency range, which is to be transmitted to the FFT analyzer 
94. In this arrangement, the FFT analyzer 94 bears a reduced burden, and assures improved analyzing accuracy. 
It will be understood from the above description that the FFT analyzer 94 and the computer 95 constitute the esti- 

45 mating means for estimating a degree of contact of the tired vehicle wheel with a road surface, and that the pre-processing 
filter 64 constitutes frequency characteristic controlling for providing an optimum frequency characteristic of the wheel 
speed signal as the wheel motion state, which is to be supplied to the estimating means. 

Referring next to Figs. 44-46, there will be described an eighth embodiment of the present invention, which includes 
the rotor 10, magnetic pickup 12, teeth 16, waveform shaper 18. computers 20 and 47, display device 66, torque detector 

so 68 and vehicle speed sensor 70, as in the first embodiment shown in Fig. 1. The functions of these elements are sub- 
stantially the same as those of the elements of the first embodiment. The present embodiment is different from the first 
embodiment in that the pre-processing filter 64 is not provided between the computers 20, 47, and that the ROM 49 of 
the computer 47 stores various control programs including a tire air pressure detecting routine as illustrated in the flow 
chart of Fig. 46. Thus, the computer 47 constitutes a disturbance observer 52. pre-processor 54, correlation calculator 

55 56. normalizer 58. basis value compensator 60, determining processor 62, and a wheel speed generator 65, as shown 
in Fig. 45. 

The rim speed calculator/compensator 45 shown in Fig. 45 is adapted to calculate the wheel speed v of each of the 
tired wheels 1 4, and adjust the calculated speed of each wheel 1 4, on the basis of the vehicle running speed (speed of 
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the vehicle body 72) detected by a suitable vehicle speed sensor 70 as indicated in Fig. 44, in the same manner as in 
the first embodiment. That is, the routines of Figs. 5 and 6 are executed by the rim speed calculator/compensator 45. 

The thus obtained wheel speed v is used by the disturbance observer 52 and the wheel speed generator 65 as 
s^own in Fig. 45. In the same manner as in the first embodiment, the disturbance observer 52 operates to estimate the 

5 disturbance w 2 = (-1/J B )T d + (AK/J B )6 RB , due to the amount of change AK of the spring constant K of the torsion 
spring 32. Similarly, the disturbance observer 52 operates to estimate the disturbance due to the change AJ B of the 
moment of inertia J B of the belt 30, and the disturbance due to the change AJ R of the moment of inertia J R of the rim 28. 

The disturbance observer 52 receives an angular velocity co R as the wheel motion state indicative of an amount of 
the motion of the tired wheel. The angular velocity co R is calculated, taking the radius R of the tire 26 into account, from 

10 the wheel speed v of the wheel 14 that is calculated and compensated by the rim speed calculator/compensator 45. 
Then, the disturbance observer 52 estimates the disturbance w 2 as indicated above, due to the amount of change AK 
of the spring constant K of the torsion spring 32, on the basis of basic values of wheel information, including a basic 
value of the moment of inertia J R of the rim 28, a basic value of the moment of inertia J B of the belt 30, and a basic value 
of the spring constant K of the torsion spring 32 provided between the rim 28 and the belt 30. Thus, the estimated 

is disturbance value w 2p is obtained, and at the same time the angular velocity « B of the belt 30 and the torsion angle e RB 
between the rim 28 and the belt 30 are estimated, to provide respective estimated values co BP e RBp . 

The pre-processor 54 operates to perform a preliminary processing operation necessary for the correlation calculator 
56 to operate. That is, the pre-processor 54 is adapted to calculate an angular acceleration a> R of the rim 28 and an 
estimated value cogp' of the angular acceleration co B \ on the basis of the detected angular velocity co R of the rim 28 and 

20 the estimated angular velocity co^ of the belt 30 obtained by the disturbance observer 52. 

The correlation calculator 56 operates to calculate correlation functions, on the basis of the estimated disturbance 
w 2p and torsion angle e RBp , by executing the flow chart of Fig. 1 1 as used in the first embodiment. Using the calculated 
correlation functions, the normalizer 58 performs a normalizing operation and calculates the amount of change AK of 
the spring constant of the torsion spring 32, from the above-indicated basic value K of the spring constant, in the same 

25 manner as in the first embodiment. As a result, the value L K = C(w 2p , 8 RBp )/C(e RBp . 6 RBp ) is obtained by the nor- 
malizer 58. 

The basic value compensator 60 operates to compensate the spring constant basic value K of the torsion spring 

32, on the basis of the above value L K stored in the L K memory of the RAM 50. The value L K is also expressed by the 

following equation: L K = (-1/J B )C 0 + AK/J B , and the amount of change AK of the spring constant is calculated accord- 
30 ing to a data table stored in the ROM 49, which represents a relationship between the value L K and the amount of change 

AK. The basic value compensator 60 is adapted to compensate the spring constant basic value K of the disturbance 

observer 52, on the basis of the calculated amount of change AK. 

When the disturbance observer 52 is operated for the first time after a power-on switch of the vehicle is turned on, 

a nominal or predetermined value is used as the basic value K of the spring constant. Namely, the disturbance observer 
35 52 uses the nominal value as the initial value of the spring constant basic value K. Once the initial value is compensated, 

however, the compensated value is used as the spring constant basic value K. The spring constant basic value K which 

has been compensated is stored in a "basic value" memory of the RAM 50. 

The compensation of the spring constant basic value K is repeated until the absolute value of the amount of change 

AK of the spring constant becomes equal to or smaller than a predetermined reference value A Kg (which is 0 or approx- 
40 imately 0, for example). Namely, the compensation is repeated until the spring constant basic value K used by the 

disturbance observer 52 becomes sufficiently close to an actual or current spring constant of the torsion spring 32. Thus, 

the basic value K obtained as a result of the compensation is estimated as the current spring constant of the torsion 

spring 32. 

In the present embodiment, the spring constant basic value K is compensated on the basis of the amount of change 
45 AK, such that the amount and direction (negative or positive) of compensation of the basic value K are determined 
depending upon the amount of change AK. However, the amount of compensation of the basic value K may be a fixed 
value irrespective of the amount of change AK, while the direction of compensation of the basic value K is determined 
depending on whether the amount of change AK is positive or negative. 

The determining processor 62 operates to read the spring constant basic value K stored in the "basic value" memory 
so of the RAM 50 as the estimated spring constant of the torsion spring 32, and estimate the air pressure P of the tire 26 
on the basis of the basic value K. When the thus estimated air pressure P is lower than a predetermined reference value 
P 0 . this means that the air pressure P of the tire 26 is abnormally low. and the display device 66 provides an indication 
informing the vehicle operator of this fact. 

While the respective operations of the disturbance observer 52. correlation calculator 56. normalizer 58, basic value 
55 compensator 60 and determining processor 62 have been described in detail, the functional relationship among these 
elements will be hereinafter described, referring to the tire air pressure detecting routine as illustrated in the flow chart 
of Fig. 46. 
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In the first step S301 , the disturbance observer 52 estimates the disturbance w 2 , the angular velocity of the belt 
30 and the torsion angle 6 RQ| on the basis of the angular velocity o R of the rim 28 received from the rim speed calcula- 
tor/compensator 45. and the driving or braking torque received from the torque detector 68. 

Step S302 is then implemented to execute the routine for calculating correlation functions for variation of the spring 

s constant K, as illustrated in the flow chart of Fig. 1 1 . More specifically, a cross-correlation function C(w 2p , e RB p) and an 
auto-correlation function C(e RB p, 8 RBp ) are respectively calculated, and the value L K is obtained by dividing the cross- 
correlation function C(w2p. 8 RBp ) by the auto-correlation function C(6 RBp , 8 RBp ). Then, the amount of change AK of the 
spring constant (an amount of change from the basic value K) is obtained on the basis of the thus obtained value 
using the data table representative of the Lk-AK relationship. 

io Step S302 is followed by step S303 to determine whether the absolute value of the amount of change AK of the 
spring constant is equal to or smaller than the reference value AK S . If an affirmative decision (YES) is obtained in this 
cycle in step S303. step S304 is then implemented to read the current spring constant basic value K from the "basic 
value" memory, as the actual spring constant. Step S305 is then implemented to estimate the current level of the air 
pressure P of the tire 26. on the basis of the spring constant read in step S304, according to a data table representative 

is of a predetermined relationship between the air pressure P and the spring constant. This data table is stored in the ROM 
49 of the computer 47. 

The control flow then goes to step S306 to determine whether the air pressure P is lower than the reference value 
Po, that is, whether the tire air pressure P is abnormally low. If a negative decision (NO) is obtained in step S306, the 
control flow goes back to step S301 . If an affirmative decision (YES) is obtained in step S306, step S307 is implemented 

20 to activate the display device 66 to provide an indication informing the vehicle operator that the air pressure P is abnor- 
mally low. In this manner, one cycle of execution of the routine of Fig. 46 is terminated. 

If the absolute value of the amount of change AK of the spring constant is larger than the reference value AKs. a 
negative decision (NO) is obtained in step S303. and step S308 is implemented to compensate the spring constant 
basic value K used by the disturbance observer 52. The compensation may be effected by adding the amount of change 

25 AK calculated in the current cycle to the basic value K, or using a predetermined relationship among the basic value K, 
amount of change AK and amount of compensation for the basic value K. After the basic value K is compensated in 
step S308. the control flow goes back to step S301. and the disturbance observer 52 estimates the above values w 2 . 
co B and 8 RB , on the basis of the compensated spring constant basic value K. If the absolute value of the amount of 
change AK is still larger than the reference value AKs, spite of the compensation effected in the previous cycle, a 

30 negative decision (NO) is obtained again in step S303. and the basic value K is compensated again in step S308. The 
compensation of step S308 is repeated until the absolute value of the amount of change AK becomes equal to or smaller 
than the reference value AKg, and an affirmative decision (YES) is obtained in step S303. 

The wheel speed generator 65 compensates the wheel speed v received from the rim speed calculator/compensator 
45, on the basis of the disturbance estimated by the disturbance observer 52, and provides an output representative of 

35 the compensated wheel speed v. 

As described above, the disturbance w 2p estimated by the disturbance observer 52 on the basis of the equation 
(1 1) is expressed by the equation (8): w 2p = (-1/J B )T d + (AK/J B )8 RB . The second term of the right member of the 
equation (8) is successively compensated by the basic value compensator 62, and therefore will not vary abruptly or at 
a high rate. This means that the second term of the right member is negligibly small as compared with the first term of 

40 the same member. Accordingly, the wheel speed generator 65 compensates the wheel speed v on the assumption that 
the estimated disturbance w^ obtained by the disturbance observer 52 according to the equation (1 1) is equal to (- 
1/J B )T d . 

Described in detail, the disturbance torque T d is obtained by multiplying the estimated disturbance w 2p = (-1/J B )T d 
by -J B , and the angular velocity Wfp of the rim 28 which is caused solely by the disturbance torque T d is estimated 
45 according to the following equation (37): 

co Rp (s) = {[D] (s[l] - [E]) 1 [F]} T d (s) (37) 

where, 

50 

[I]: unit (identity) matrix, 

s : Laplace operator, 

«>Rp(s): value obtained by Laplace transform from the estimated angular velocity copp, 

T d (s): value obtained by Laplace transform from the disturbance torque T d . 
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The [D], [E] and [F] are vectors and a matrix expressed by the following equations: 
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The estimated angular velocity or p is one of the components causing the speed variation of the tired wheel 14, 
which one component is caused by the disturbance given to the wheel 1 4 from the road surface irregularity. The estimated 
angular velocity a>Rp is converted into the corresponding peripheral speed of the wheel 14, and the wheel speed v 
received from the rim speed calculator/compensator 45 is compensated by the amount corresponding to the calculated 
peripheral speed, whereby the noise due to the disturbance from the road surface is eliminated from the wheel speed v. 

It will be understood from the foregoing description of the present embodiment that the rotor 10, magnetic pickup 
12, waveform shaper 18 and the calculator of the rim speed calculator/compensator 45 for calculating the wheel speed 
v cooperate to provide a detecting device for detecting the wheel motion state. It will also be understood that a portion 
of the computer 47 assigned to execute step S301 of the routine of Fig. 46 constitutes a disturbance observer, and a 
portion of the computer 47 assigned to execute step S302 of the same routine constitutes means for estimating an 
amount of change from the basic value to the current value of the wheel information. Further, portions of the computer 
47 assigned to execute steps S303 and S308 of the same routine constitute basic value compensating means, while 
portions thereof assigned to execute steps S304 and S305 of the same routine constitute wheel information estimating 
means 

The tired wheel 14 is dynamically simulated as a system or model as illustrated in Fig. 47, taking account of a 
damper between the rim 28 and the belt 30. Namely, the dynamic model is designed such that the mutually rotatable 
rim 28 and the belt 30 are connected to each other, by the torsion spring 32 and a damper 130 which are disposed in 
parallel with each other. 

The model of Fig. 47 may be simplified into a so-called "two-inertial model" as illustrated in Fig. 48. 
The reference characters used in Fig. 47 represent the following: 



m R : equivalent inertial mass (weight) of the rim 28, 

m B : equivalent inertial mass (weight) of the belt 30, 

k w : spring constant of the torsion spring 32, 

D w : damping coefficient of the damper 1 30, 

x R : equivalent linear displacement of the rim 28, 

x B : equivalent linear displacement of the belt 30, 

x RB : equivalent relative linear displacement of the rim 28 and belt 30, 

F d : equivalent disturbance received by the tire 26 from the road surface. 

The equivalent inertial masses (weights) m R , m B of the rim 28 and belt 30 in the model of Fig. 48 equivalents 
correspond to the moments of inertia J R , J B in the model of Fig. 47, respectively, while the equivalent linear displacements 
x R , x B of the rim 28 and belt 30 in the model of Fig. 48 equivalents correspond to the integrals of the angular velocities 
co R , co B in the model of Fig. 47, respectively. The equivalent relative linear displacement x RB in the model of Fig. 48 
equivalents corresponds to the torsion angle 9 RB in the model of Fig. 47. Further, the equivalent disturbance F d in the 
model of Fig. 48 equivalent^ corresponds to the disturbance torque T d in the model of Fig. 47. There will be described 
a ninth embodiment of the present invention according to the model of Fig. 48. 
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The dynamic model of Fig. 48 is represented by the following state equation (38): 




75 The above equation (38) does not include a parameter corresponding to the driving or braking torque T-i which 

acts on the rim 28 based on an engine output or a braking force, because the equation (38) is formulated to include only 
variable parameters, with a particular attention drawn to vibrating movements of the model of Fig. 48. In this respect, 
the driving or braking torque T 1 can be considered a fixed value as distinguished from the variable parameters, such as 
the equivalent linear displacements x R , x B , first derivatives x R \ x B * and second derivatives x R ", x B " of those values x R> 

20 x B , spring constant K w , and damping coefficient D w 

If the spring constant K w and damping coefficient D w are changed to (K w + AK W ) and (D w + AD W ), respectively, 
due to a change of the air pressure P of the tire 26, the state equation (38) can be converted into the following equation 
(39), which represents an equivalent state of the model in which a disturbance acts on the model placed in its normal 
state: 
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If the equivalent disturbance F d is also treated as an unknown disturbance value that cannot be measured, the 
disturbance to be estimated is represented by the following equation (40): 
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If the disturbance represented by the second elements of the right member of the above equation (40) is to be 
estimated, the disturbance w 2 is represented by the following equation (41): 



w 2 = (AD w /m B )(x „' - x B ) + (AK w /m B )x RB + F d /m B + n 



(41) 



10 



where, "n" represents a term representative of an error which arises because only the second elements are used 
to estimate the disturbance w 2 . 

Therefore, the following state equation (42) represents the model when the air pressure of the tire 26 is changed: 

fx* " 1 r-Dw /m R Dw /m R -Kw /m R x R ' 1 

X e " = Dw /ibb "Dw /m B Kw Mb X b ' 

LXrb' J L l - l o Xrb J 
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The wheel .nformat.on estimating apparatus according to the present embodiment is based on the above analysis 
and includes a disturbance observer 160 to estimate the disturbance w 2 for obtaining the amount of change of the air 
pressure P of the tire 26. The present wheel information estimating apparatus is constructed as described below 

The apparatus includes a displacement velocity detecting device 1 40 as illustrated in the functional block diaphragm 
of Fig. 49. The device 140 is adapted to detect an equivalent linear displacement velocity x R ' of the rim 28 and is 
connected to a computer 150. Thedevice 1 40 operates to determine the angular velocity <o R of the rim 28 by magnetically 
detecting the passage of the teeth 1 6 formed along the outer periphery of the rotor 1 0. The equivalent linear displacement 
velocity x R - is calculated on the basis of the determined angular velocity co R . As shown in Fig. 50. the computer 150 
incorporates a CPU 152 as a central processing unit, a ROM 154 as a first memory device, and a RAM 156 as a second 
memory device. The computer 1 50 provides the above-indicated disturbance observer 160, a constant change calculator 
162. a basic value compensator 164 and a determining processor 166, as indicated in Fig. 49 

The disturbance observer 160 is identical in construction with the disturbance observer 52 used in the previous 
embodiments, except for the following: 

If the disturbance to be estimated is approximated as w 2 ' = 0, the expanded system of the linear system represented 
by the above equation (42) is represented by the following equation (43) : 
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(43) 



Of the physical values in the above equation (43), only the equivalent linear displacement velocity x R of the rim 
28 can be detected. 
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The disturbance observer 1 60 can be formulated in the same manner as described above, by re-defining the matrices 
in the above equation (43) as follows: 
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20 In essence, the present embodiment is adapted such that the disturbance observer 160 estimates the equivalent 

linear displacement velocity x B and the equivalent relative linear displacement x RB of the linear system relating to the 
rotation of the tired wheel, on the basis of at least the equivalent linear displacement velocity x R ' of the rim 28 detected 
by the displacement velocity detector 140, equivalent inertia mass m R of the rim 28, equivalent inertia mass m B of the 
belt 30, spring constant basic value K w of the torsion spring 32 between the rim 28 and the belt 30. and damping 

25 coefficient basic value K w of the damper 130 between the rim 28 and the belt 30. The disturbance observer 160 also 
estimates the disturbance w 2 in the linear system. It will be understood that the above-indicated equivalent linear dis- 
placement velocity x R ' is the wheel motion state indicative of the motion of the tired wheel 1 4. and the equivalent inertia 
masses m R , m B , spring constant basic value Kw, and damper coefficient basic value Dw are basic values indicative of 
wheel information relating to various conditions of the wheel 14. 

30 The disturbance w 2 thus estimated by the disturbance observer 160 is applied to the constant change calculator 
1 62, which is adapted to calculate an amount of change AK W of the spring constant K w and an amount of change ADw 
of the damping coefficient D w . on the basis of the disturbance w 2 . 

As described above. w 2 is represented by the following equation: 
w 2 = (AD w /m B )(x R '-x B ') + (AK w /m B )x RB + F d /m B + n . This equation is used to calculate the amount of change 

35 AK W and the amount of change AD W by the least square method. Explained more specifically, the amounts of change 
AK W and AD W are determined such that a sum of the squares expressed by the following equation (44) is minimized, 
that is, such that values obtained by partial differentiation of the sum S by the respective amounts of change AK W and 
• AD W are both zeroed. 
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The amounts of change AK W and AD W are calculated using the following equation (45): 
A Kw 
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where, 

R 1(i) " x RB(i) 
R 2(i) = X R (j) ' 



x b'o) 



The thus obtained amounts of change AK W and AD W are applied to the basic value compensator 164, which 
compensates or adjusts the spring constant basic value K w and the damping coefficient basic value D w , according to 
the amounts of changes AK W and AD W , respectively. As in the previous embodiment, the compensation is repeated until 
the absolute value of the amount of change AK W of the spring constant becomes equal to or smaller than a predetermined 
reference value AK WS . The compensated basic values K w and D w are stored in a "basic value" memory of the RAM 1 56. 

The compensated basic values K w and D w of the spring constant and damping coefficient are applied to the deter- 
mining processor 166, assuming that the values KW and DW represent actual spring constant and damping coefficient. 
The determining processor 1 66 estimates the present or current value of the air pressure P, on the basis of the present 
values of the spring constant and the damping coefficient, according to a predetermined relationship between the spring 
constant and damping coefficient, and the air pressure P of the tire 26. This relationship is represented by a data table 
stored in the ROM 154 of the computer 150. Finally, the processor 166 determines whether the estimated air pressure 
P is lower than a predetermined reference value Po. and activates a display device 1 68 to provide an indication that the 
air pressure P of the tire is abnormal, if the estimated air pressure P is lower than the reference value Po- 

The operation described above is performed according to a tire air pressure detecting routine stored in the ROM 
154 of the computer 150. This routine will be described in detail by reference to the flow chart of Fig. 51 . 

The routine is initiated with step S31 1 to set the integer V to "1 ", and reset the amount of change AK W(l) and AD w(i) 
to "0". Then, step S312 is implemented to read the equivalent linear displacement velocity x R (i) of the rim 28 from the 
displacement velocity detector 140, and store the velocity Xr'^ in the RAM 156. 

Step S31 2 is followed by step S31 3 to calculate the equivalent linear displacement acceleration Xr"^ of the rim 28, 
which is a difference between the present velocity value x R '(j) and the last velocity value x R '( M ) previously stored in the 
RAM 156. The calculated equivalent linear displacement acceleration x R " (j) is also stored in the RAM 156. 

The control flow then goes to step S314 in which the disturbance observer 160 estimates the disturbance w 2 (j), and 
the equivalent linear displacement velocity x B (j) and equivalent relative linear displacement x RB(i ) of the belt 30, on the 
basis of the equivalent linear displacement velocity xr'^ and the equivalent linear displacement acceleration x R " (j j. The 
estimated values w 2 (j), x B '(j) and x RB (|) are stored in the RAM 156. 

While the equivalent linear displacement acceleration x R "(j) as well as the equivalent linear displacement velocity 
x R ' (i) is used to estimate the disturbance w 2 {\) and others, only the equivalent linear displacement velocity x R ' (i) may be 
used by the disturbance observer 160. 

The disturbance observer 160 receives the equivalent linear displacement velocity and acceleration x R '(j), x R "(j) only 
after these inputs x R ' w , x R " (i ) are filtered by a high-pass filter (not shown), so that only the variable components of the 
velocity and acceleration x R ' ( j), x R ^ are used by the observer 160. 

Then, step S315 is implemented to determine whether the integer V is equal to or larger than a predetermined 
value "N" or not. If not, a negative decision (NO) is obtained in step S315, and step S316 is implemented to increment 
the integer T to repeat the above steps S312-S315. With these steps S312-S316 repeated the predetermined number 
"N" of times, a total of "N" sets of xr'^, x B ' (i) , w 2 (j) and x RB(j ) are stored in the RAM 156, and an affirmative decision 
(YES) is obtained in step S315. In this case, step S317 is implemented to calculate the amounts of change AK w(i) and 
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AD W(j) . More specifically, these amounts of change AK W(i) and AD W(i) are calculated by the least square method on the 
basis of the M N H sets of data indicated above. Step S317 is followed by step S318 to determine whether the absolute 
value of the amount of change AK w(j) is equal to or smaller than the reference value AK wS or not. If an affirmative 
decision (YES) is obtained in step S318, the control flow goes to step S319. 

5 In step S31 9, the current level of the air pressure P is estimated, on the basis of the spring constant basic value K w 

and the damping coefficient basic value Dw. which are retrieved from the "basic value" memory of the RAM 156, and 
which are considered to represent the actual spring constant and damping coefficient. Step S319 is followed by step 
S320 to determine whether the air pressure P is lower than the reference value Po or not. If not, the display device 168 
is not activated to inform the vehicle operator of the abnormal air pressure P. In this case, a negative decision (NO) is 

10 obtained in step S320. and step S322 is then implemented to reset the integer T to 1 . Then, the control flow goes back 
to step S312, and subsequent steps are executed in the next cycle of execution of this routine. 

If the air pressure P is lower than the reference value Po. an affirmative decision (YES) is obtained in step S320, 
and step S321 is implemented to activate the display device 168 to provide an indication to inform the vehicle operator 
that the air pressure P is abnormal. In this case, the execution of the routine is terminated with step S321 . 

75 If the absolute value of the present amount of change AK W(i) of the spring constant is larger than the reference value 
AK W s. a negative decision (NO) is obtained in step S318, and step S323 is implemented to compensate the spring 
constant basic value K w and the damping coefficient basic value D w used by the disturbance observer 160, on the basis 
of the currently obtained amounts of changes AK W (j) and AD W(i ). This step S323 is equivalent to step S308 of the routine 
of Fig. 46 of the previous embodiment. After executing step S323, the control flow goes back to step S312, and steps 

20 S312-323 are repeated until an affirmative decision (YES) is obtained in step S318, that is, until the absolute value of 
the amount of change AK W(I) becomes equal to or smaller than the reference value AK WS . 

While only the amount of change AK W ^ is compared with the reference value AK W s in step S318, the amount of 
change AD W(i ) as well as the amount of change AK W (q may be compared with the reference value ADws For example, 
the compensation of the basic values K w , Dyy may be repeated until the absolute values of the amounts of change AK W(i) 

25 and AD W(i) become equal to or smaller than the respective reference values AK WS and AD WS , and an affirmative decision 
(YES) is obtained in step S3 18. 

While the "basic value" is a predetermined nominal value in the present embodiment, the basic value according to 
the present invention may be 0 or any value other than the nominal value. 

It will be understood from the above description of the present embodiment that the rotor 1 0. magnetic pickup 12, 

30 waveform shaper 1 8. the calculator of the rim speed calculator/compensator 45 for calculating the wheel speed cooperate 
to provide a detecting device for detecting the wheel motion state as indicated above. It will also be understood that 
portions of the computer 47 assigned to execute steps S31 2 and S31 4 of the routine of Fig. 51 constitute a disturbance 
observer, and that a portion of the computer 47 assigned to execute step S317 of the same routine constitutes means 
for estimating an amount of change from the basic value to the current value of the wheel information as indicated above. 

35 Further, portions of the computer 47 assigned to execute steps S323 and S318 of the same routine constitute basic 
value compensating means, while a portion thereof assigned to execute step S3 1 9 of the same routine constitutes wheel 
information estimating means. 

It will be understood that the tire air pressure detecting routine of Fig. 51 is formulated such that the amounts of 
change AK W and AD^ are calculated each time the predetermined "N" set of data x R ' (i) , x B ' (j) , w 2( i) and x RB(i) have been 

40 stored in the RAM 156. However, the routine may be modified such that the amounts of change AK W and AD W are 
calculated each time a new set of those data has been stored in the RAM 156 after "N" sets of the data have been 
stored. In this case, the above equation (45) is modified so as to update the "N" sets of data such that the first stored 
set of data (oldest set of data) is replaced by the last stored set of data (newly obtained set of data). An example of such 
modified routine is illustrated in the flow chart of Fig. 52. This modified routine will be described, with a brief explanation 

45 of the steps similar to those of the routine of Fig. 51 . 

The routine of Fig. 52 is initiated with step S341 to set the integer V to "1 and zero not only the amounts of change 
AK W (j) and AD w( q but also matrices Lq and L^. The matrix Lq is represented by the left term of the right member of the 
above equation (45), while the matrix Lm is represented by the right term of the right member of the equation (45). 
Step S341 is followed by step S342 to calculate the equivalent linear displacement velocity x R ' (j j on the basis of the 

so output of the displacement velocity detector 140. Step S343 is then implemented to calculate the equivalent linear dis- 
placement acceleration x R " (i) . Step S343 is followed by step S344 in which the disturbance observer 160 estimates the 
disturbance w 2 ^, the equivalent linear displacement velocity xq'q and the equivalent relative linear displacement x RB(i ), 
on the basis of the variable components of the velocity and acceleration x R ' (i)l x R " ( j). 

The control flow then goes to step S345 to add to the present value of the matrix Lo an increment based on the 

55 present values of the equivalent linear displacement velocities xr^). x B(i) and the equivalent relative linear displacement 
x RB(i) Step S346 is then implemented to add to the present value of the matrix L N an increment based on the present 
values of x R ' (i) , x B ' ( i), x RB(i j and the present value of the disturbance w 2 (j). Thus, steps S345 and 346 are provided to 
update the matrices L D and Ln. Then, the control flow goes to step S347 to calculate the amounts of change AK W (j) and 
AD W(j ) by inserting Lq and L N in the equation (45). Step S347 is followed by step S348 to compensate the spring constant 
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basic value K w and the damping coefficient basic value D w used by the disturbance observer 160, on the basis of the 
calculated amounts of change AK W(i) and AD W(i) . The content of this step S348 is the same as that of steps S3 18 and 
S323 of the routine of Fig. 51 . 

The control flow then goes to step S349 to increment the integer V, and then to step S350 to determine whether 

5 the air pressure P is abnormal or not. This determination is made on the assumption that the present basic values K w 
and Dw represent actual spring constant and damping coefficient, respectively. If an affirmative decision (YES) is 
obtained in step S350, step S351 is implemented to activate the display device 66 to inform the vehicle operator of the 
abnormal air pressure P. The control flow then goes back to step S342. 

While the present invention has been described in detail above in its presently preferred embodiments, it is to be 

10 understood that the present invention is not limited to the details of the illustrated embodiments, but may be modified 
with changes, modifications and improvements, which may occur to those skilled in the art, without departing from the 
scope of the invention defined in the following claims. 

A wheel information estimating apparatus, which includes: a detecting device (10, 12, 18, 56, 140) for detecting a 
wheel motion state relating to a movement of a tired wheel of a motor vehicle; an estimating device (52, 53, 82, 84, 86, 

15 94, 95, 160, 162) for estimating wheel information relating to the tired wheel, on the basis of the wheel motion state 
detected by the detecting device; and a modifying device (64, 80, 92, 60, 164) for modifying an element relating to 
estimation of the wheel information, so as to improve an accuracy of estimating the wheel information. The element to 
be modified is selected from a movement of a center of rotation of the tired wheel as viewed in a plane in which the 
wheel is rotated, a signal applied from the detecting device to the estimating means, and an internal parameter used by 

20 the estimating means for estimating the wheel information. 

Claims 

1. A wheel information estimating apparatus, including a detecting device (10, 12, 18. 45, 140) for detecting a wheel 
25 motion state relating to a movement of a tired wheel (1 4) of a motor vehicle, and estimating means (52, 53, 82, 84, 

86. 94. 95, 160. 162) for estimating wheel information relating to the tired wheel, on the basis of said wheel motion 
state detected by said detecting device, characterized in that: 

modifying means (64, 80, 92, 60, 164) is provided for modifying an element relating to estimation of said 
wheel information, so as to improve accuracy of estimating the wheel information, said element being at least one 
30 of: a movement of a center of rotation of the tired wheel as viewed in a plane in which the wheel is rotated; a signal 
applied from said detecting device to said estimating means; and an internal parameter used by said estimating 
means for estimating the wheel information. 

2. A wheel information estimating apparatus according to claim 1 , wherein said modifying means comprises frequency 
35 characteristic modifying means (64, 80, 92, 94, 95) for providing an optimum frequency characteristic of said wheel 

motion state to be supplied to said estimating means, by increasing a ratio of a strength of a selected one of a 
plurality of frequency components of said wheel motion state, said selected one being within a selected frequency 
range, to a strength of the other frequency components outside said selected frequency range. 

40 3. A wheel information estimating apparatus according to claim 2, wherein said modifying means modifies said signal 
applied from said detecting device to said estimating means, said signal indicative of said wheel motion state, said 
frequency characteristic modifying means comprising a pre-processing filter (64, 80, 92) disposed between said 
detecting device and said estimating means, for selecting a component of said signal received from said detecting 
device, which component is within said selected frequency range, and supplying the selected component to said 

45 estimating means. 

4. A wheel information estimating apparatus according to claim 2 or 3, wherein said frequency characteristic modifying 
means comprises means (90, S104, S105, S1 14, S1 15, S1 16) for changing said selected frequency range, on the 
basis of at least one of an estimated value of said wheel information obtained by said estimating means, and a 

so related value of said estimated value. 

5. A wheel information estimating apparatus according to any one of claims 2-4, wherein said frequency characteristic 
modifying means comprises means (90, S126, 127) for changing said selected frequency range, so that an amount 
of variation of an estimated value of said wheel information obtained by said estimating means is not greater than 

55 a predetermined value. 

6. A wheel information estimating apparatus according to any one of claims 2-5, wherein said frequency characteristic 
modifying means comprises a geometry of a suspension which connects the tired wheel and a vehicle body such 
that the wheel and the vehicle body are displaceable relative to each other, said geometry being designed such that 
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a center of the tired wheel provides an optimum locus as seen in a transverse direction of the motor vehicle, so as 
to improve accuracy of estimating the wheel information, said optimum locus of the center of the tired wheel being 
drawn upon a relative movement between the wheel and the vehicle body. 

7. A wheel information estimating apparatus according to claim 6, wherein said suspension comprises an upper arm 
(104) and a lower arm (106) which extend in a longitudinal direction of the motor vehicle, and cooperate with each 
other to connect said tired wheel to said vehicle body, such that the wheel and the vehicle body are displaceable 
relative to each other, said geometry of the suspension being designed such that said upper arm extends substan- 
tially in parallel with said lower arms as seen in the transverse direction of the vehicle, so that said locus of the center 
of the tired wheel is drawn substantially in a vertical plane which extends in the transverse direction of the motor 
vehicle. 

8. A wheel information estimating apparatus according to any one of claims 2-7, wherein said estimating means com- 
prises (a) a disturbance observer (52) for estimating a disturbance acting on the tired wheel, on the basis of a basic 
value of said wheel information, and said wheel motion state received from said detecting device, and (b) means 
(56, 58, 62) for estimating an amount of change from said basic value to a value indicative of an actual state of said 
wheel information, on the basis of said disturbance estimated by said disturbance observer. 

9. A wheel information estimating apparatus according to any one of claims 2-7, wherein said estimating means com- 
prises means for estimating said wheel information, on the basis of a frequency of a selected one of a plurality of 
frequency components of said wheel motion state received from said detecting device, said selected one of the 
frequency components having substantially the greatest strength within a selected frequency range. 

10. A wheel information estimating apparatus according to any one of claims 2-7, wherein said detecting device com- 
prises a wheel speed detector (10, 12. 18, 45) for detecting, as said wheel motion state, a wheel speed which is a 
rotating speed of the tired wheel, and supplying said wheel speed to said estimating means, said estimating means 
estimating a degree of contact of the tired wheel with a road surface as said wheel information, on the basis of a 
strength of a selected one of a plurality of frequency components of said wheel speed received from said wheel 
speed detector, said selected one of said frequency components being within a selected frequency range. 

11. A wheel information estimating apparatus according to any one of claims 2-10, wherein said estimating means 
obtains a plurality of provisional estimated values of at least one of said wheel information and a parameter relating 
to the wheel information, so as to obtain a current final estimated value of said wheel information, on the basis of 
said plurality of provisional estimated values, a number of said provisional estimated values being used for at least 
one of re-estimation of said current final estimated value and estimation of a next final estimated value of said wheel 
information, said number of the provisional estimated values being variable depending upon whether an amount of 
change from a previously estimated value of said wheel information to a currently estimated value thereof is larger 
than a predetermined value. 

1 2. A wheel information estimating apparatus according to claim 1 , wherein said modifying means modifies said internal 
parameter used by said estimating means, said estimating means comprising (a) a disturbance observer (52) for 
estimating a disturbance acting on the tired wheel, on the basis of at least said wheel motion state detected by said 
detecting device, and a basic value of said wheel information, said disturbance observer using said basic value as 
said internal parameter, (b) change amount estimating means (56, 58, 62) for estimating an amount of change from 
said basic value to a value indicative of an actual state of said wheel information, on the basis of said disturbance 
estimated by said disturbance observer, (c) basic value compensating means (60) as said modifying means for 
compensating said basic value of the wheel information to provide a compensated wheel information basic value, 
on the basis of said amount of change estimated by said change amount estimating means, so that the estimated 
amount of change is not larger than a predetermined value; and (d) wheel information estimating means for esti- 
mating said compensated wheel information basic value as said wheel information. 

13. A wheel information estimating apparatus according to claim 12. wherein said disturbance observer assumes a 
model of the tired wheel having a rim (28) and a belt (30) which are connected to each other by a torsion spring 
(32). such that said rim and said belt are rotatable relative to each other, said disturbance observer (52, S301) 
estimating an angular velocity a> B of the belt, a torsion angle 6 RB between the rim and the belt, and a disturbance 
w 2 , on the basis of a detected value of an angular velocity o> R of the rim as said wheel motion state, and basic values 
of a moment of inertia J R of the rim, a moment of inertia J B of the belt, and a spring constant K of said torsion spring, 
as said wheel information basic value, said change amount estimating means (S302) estimating an amount of 
change AK from the basic value of said spring constant K to a value indicative of an actual spring constant, on the 
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basis of estimated values of said angular velocity <ob of the belt, said torsion angle e RB , and said disturbance w 2 , 
said basic value compensating means (S303, S308) compensating the basic value of said spring constant K, so 
that the estimated amount of change AK is not larger than a predetermined value, said wheel information estimating 
means (S304, S305) estimating the basic value of the spring constant K as said wheel information, when said 
amount of change AK is not larger than said predetermined value. 

14. A wheel information estimating apparatus according to claim 13, wherein said disturbance observer defines said 
disturbance w 2 as represented by the following equation: 

w 2 =(-1/J B )T d + (AK/J B )8 RB 

where, 

T d : disturbance torque applied from a road surface to said belt of the tired wheel, 

and wherein said change amount estimating means estimates said amount of change AK according to a relationship 
between the estimated values of said disturbance w 2 and said torsion angle 6 RB . 

15. A wheel information estimating apparatus according to claim 14, wherein said change amount estimating means 
calculates a normalized value obtained by dividing a cross-correlation function of the estimated values of said dis- 
turbance w 2 and said torsion angle e RB , by an auto-correlation function of the estimated value of said torsion angle 
e RB , and estimates said amount of change AK, on the basis of said normalized value as a variable representative 
of a relationship between said disturbance w 2 and said torsion angle 9 RB . 

16. A wheel information estimating apparatus according to claim 12, wherein said disturbance observer assumes an 
equivalent model of the tired wheel which is equivalent to a basic model in which a rim (28) is coupled with a belt 
(30) through a torsion spring (32) and a damper (130) which are connected in parallel with each other, such that 
said rim and said belt are rotatable relative to each other, said equivalent model being designed such that a rim is 
coupled with a belt through a torsion spring and a damper which are connected in parallel with each other, such 
that said rim and said belt are displaceable relative to each other in one linear direction, said disturbance observer 
(160, S312, S314) estimating an equivalent linear displacement velocity x B ' of the belt, an equivalent relative linear 
displacement x RB between the rim and the belt, and a disturbance w 2 , on the basis of a detected value of an equiv- 
alent linear displacement velocity x R ' of said rim, as said wheel motion state, and basic values of an equivalent 
inertial mass m R of the rim, an equivalent i nertial mass m B of the belt, a spring constant K w and a damping coefficient 
D w , as said wheel information basic value, said change amount estimating means (S317) estimating an amount of 
change AK W from the basic value of said spring constant K w to a value indicative of an actual spring constant, and 
an amount of change AD W from the basic value of said damping coefficient D w to a value indicative of an actual 
damping coefficient, on the basis of estimated values of said equivalent linear displacement velocity x B of the belt, 
said equivalent relative linear displacement x RB and said disturbance w 2 , said basic value compensating means 
(S323, S31 8) compensating the basic values of the spring constant K w and the damping coefficient D w , so that the 
estimated amounts of change AK W and AD W are not larger than respective predetermined values, said wheel infor- 
mation estimating means (S31 9) estimating the basic values of the spring constant K w and the damping coefficient 
D w as the actual spring constant and damping coefficient, when said estimated amounts of change AK W and AD W 
are not larger than said respective predetermined values. 

17. A wheel information estimating apparatus according to claim 16, wherein said disturbance observer defines said 
disturbance w 2 as represented by the following equation: 

w 2 = (AD w /m B )(x R ' - x B ') + (AK w /m B )x RB + F d /m B 

where, 

F d : disturbance received by said belt from a road surface, 

and wherein said change amount estimating means estimates the amounts of change AD W and AK W on the basis 
of estimated values of said disturbance w 2 , said equivalent linear displacement velocity x B of the belt, and said 
equivalent relative linear displacement x qB , and said detected value of the equivalent linear displacement velocity 
x R ' of the rim. 
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18. A wheel information estimating apparatus according to claim 17, wherein said wheel information estimating means 
estimates said amount of change ADw and said amount of change AK W . on the basis of the estimated values of the 
disturbance w 2 , the equivalent linear displacement velocity x B ' of the belt, and the equivalent relative linear displace- 
ment x RB , and the detected value of the equivalent linear displacement velocity x R * of the rim, such that respective 
estimated values of said amount of change AD W and said amount of change AK W are determined so as to substan- 
tially minimize a sum of the squares of a difference between the estimated value of the disturbance w 2 and a the- 
oretical approximate value thereof. 

19. A wheel information estimating apparatus according to any one of claims 12, 13 and 16, further comprising tire 
pressure estimating means for estimating an air pressure P of a tire attached to said tired wheel, on the basis of 
said wheel information estimated by said wheel information estimating means. 

20. A wheel information estimating apparatus according to claim 1 , wherein said modifying means (64, 80, 92, 60) 
performs at least one of two functions, one of which is to convert a first frequency characteristic of said wheel motion 
state detected by said detecting device, into a second frequency characteristic which improves accuracy of estimat- 
ing said wheel information, and apply said second frequency characteristic to said estimating means, the other of 
said two functions being to convert said internal parameter used by said estimating means, into a modified parameter 
which improves accuracy of estimating said wheel information, during operation of the estimating means. 
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